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ABSTRACT 


A  linearized  mathematical  model  is  developed  which  simu- 
lates the  dynamic  response  of  the  Navy  F/A-18  for  the  study 
of  flight  control  reconfiguration.   The  aircraft  is  modeled 
as  a  multi-input  multi-output,  sampled  data,  closed  system, 
which  couples  the  dynamics  of  the  flight  control  system  to 
the  aircraft  linearized  small  perturbation  equations.   The 
discrete  time,  state  variable  equations  for  the  system  are 
then  formulated.   A  computer  program  is  developed  which  will 
compose  the  model  matricies  and  compute  the  response  of  the 
aircraft  to  stick  and  rudder  inputs. 

To  study  flight  control  reconfiguration,  the  model  allows 
individual  actuation  of  either  a  left  or  right  control  sur- 
face.  Aircraft  response  to  the  actuation  loss  of  either  the 
left  or  right  stabilator  is  simulated  in  the  program.   The 
program  is  designed  to  implement  the  r econ f igurable  control 
mixer,  currently  under  study  for  the  Sel f -Repai r ing  Digital 
Flight  Control  System. 

The  computer  simulation  was  written  in  VS  FORTRAN.   A 
copy  of  the  program  and  simulation  results  are  included  in 
the  a  ppendices . 
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I  .   INTRODUCTION 

Tactical  aircraft  face  airborne  and  ground  based  threats 
which  continue  to  grow  in  number  and  capability.   In  future 
combat  engagements  the  lethality  of  these  systems  will  pose 
a  considerable  threat  to  aircraft  survivability.   In  addition 
to  combat  losses,  the  loss  of  aircraft  in  battle  damage  re- 
pair, or  awaiting  repair,  will  significantly  reduce  our  tac- 
tical forces.   Projections  on  the  survivability  of  NATO 
forces  [Ref.  1]  during  the  initial  days  of  engagement  indi- 
cate 68%  of  the  tactical  aircraft  will  be  out  of  action 
after  the  third  day  of  battle.   Of  this  total,  22%  will  be 
lost  in  combat  and  46%  lost  while  in  repair  or  awaiting  re- 
pair.  This  is  indicated  in  Fig.  1.1  which  comes  from  Refer- 
ence 1.   It  is  clear  from  projections  such   as  this,  if  our 
forces  are  to  remain  a  superior  threat  to  the  enemy,  contin- 
ued emphasis  must  be  given  to  reduce  the  combat  vulnerability, 
and  increase  the  reliability  and  maintainability  of  our  tac- 
tical aircraft. 

To  improve  the  combat  effectiveness  of  future  tactical 
aircraft,  the  United  States  Air  Force  initiated  the  Self 
Repairing  Flight  Control  System,  Reliability  and  Maintain- 
ability Program.   Reference  2  outlines  the  program  plan  and 
goals.   An  Air  Force  sponsored  study  [Ref.  1]  showed  that 
significant  improvements  in  aircraft  survivability  and 
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reliability  can  be  achieved  with  a  self  repairing  digital 
flight  control  system. 

One  of  the  principle  features  of  a  self  repairing  system 
is  flight  control  reconfiguration.   A  r eco nf i gu r ab le  flight 
control  system  automatically  counters  aircraft  loss  of  con- 
trol due  to  impairment,  or  loss,  of  a  control  effector.   The 
fundamental  idea  is  to  utilize  the  existing  undamaged  effec- 
tors to  preserve  the  normal  flying  qualities  of  the  unimpaired 
aircraft.   The  approach  outlined  in  Reference  1  is  to  design 
a  reconfigurable  control  mixer  to  be  placed  between  the  flight 
control  laws  and  the  control  effectors.   In  the  face  of  air- 
craft impairment  the  mixer  would  reallocate  control  commands 
to  the  unimpaired  effectors  so  that  flight  critical  pitch, 
yaw,  and  rolling  moments  would  be  preserved  to  the  maximum 
extent  possible.   Using  the  control  mixer  concept,  the  exist- 
ing flight  control  laws  would  not  be  altered. 

Although  the  development  of  the  self  repairing  system  is 
intended  for  the  advanced  tactical  aircraft,  it  is  possible 
that  the  mixer  could  be  implemented  in  existing  airframes. 
This  could  be  done  with  a  control  reconfiguration  module 
interposed  between  the  aircraft  flight  control  computers  and 
fly-by-wire  actuators  [Ref.  1]. 

For  this  thesis  the  McDonnell  F/A-18  was  chosen  to  study 
the  reconfigurable  control  law  concept  utilizing  the  control 
mixer  as  described  above.   This  choice  was  made  based  on  the 
following  assumptions  given  in  Reference  1  for  implementation 
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of  the  reconfiguration  control  law  process  in  an  existing 
ai  r frame  : 

1)  "The  effector  complement  provides  redundant  effector 
systems  and  surplus  control  power  for  each  flight 
critical  control  force  and  moment." 

2)  "The  flight  control  is  a  full  authority,  fly-by- wire 
digital  flight  control  system." 

3)  "A  control  law  design  exists  which  has  been  carried 
out  for  the  unimpaired  airplane,  but  which  is  suf- 
ficiently robust  that  only  first  order  impairment- 
induced  stability  derivative  changes  need  be 
accounted  for  in  a  drop-in  reconfigurable  mixer." 

Based  upon  these  assumptions  it  is  felt  that  the  F/A-18  is 

well  suited  for  the  reconfigurable  control  law  study. 

In  this  thesis  a  linearized  mathematical  model  was  de- 
veloped which  simulates  the  dynamic  response  of  the  F/A-18 
aircraft  to  stick  and  rudder  inputs.   The  model  includes  the 
control  laws  for  the  longitudinal,  lateral,  and  directional 
axes  for  the  cruise  phase  of  flight.   The  fly-by-wire  actua- 
tors and  sensor  dynamics  were  modeled  and  joined  to  the  air- 
frame linearized,  small  perturbation  model.   The  perturbation 
model  was  obtained  from  the  Simulation  Control  and  Technology 
Group,  Flight  Systems  Branch,  Strike  Aircraft  Test  Direc- 
torate, at  the  Naval  Air  Test  Center.   A  computer  program 
was  developed  which,  with  existing  software  at  the  Naval 
Postgraduate  School,  will  simulate  F/A-18  longitudinal,  lat- 
eral, and  directional  response  to  stick  and  rudder  inputs. 

The  model  was  designed  to  implement  the  control  mixer 
gain  concept  outlined  above.   Future  work  at  NPS  will  utilize 
the  simulation  program  to  develop  algorithms  for  determining 
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the  control  mixer  gain  matrix.   In  addition  the  program  will 
be  used  to  study  modern  control  augmentation  systenrs   and  air- 
craft stability  and  control. 
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II .   MODEL  METHODOLOGY 

A.   INTRODUCTION 

This  chapter  discusses  the  methods  and  assumptions  used 
to  formulate  a  mathematical  model  which  simulates  the  dynam- 
ic response  of  the  F/A-18  aircraft.   The  flight  control  sys- 
tem is  described,  including  the  simplifying  assumptions  used 
to  develop  the  control  system  model.   An  overview  of  the  com- 
plete system,  which  couples  the  control  system  model  with 
the  airframe  small  perturbation  model,  is  then  given  with  a 
brief  description  of  each  functional  component. 


B.   FLIGHT  CONTROL  SYSTEM  DESCRIPTION 

A  detailed  description  of  the  F/A-18  control  system  and 
theory  of  operation  can  be  found  in  the  Flight  Control  Sys- 
tem Design  Report  by  McDonnell  Aircraft  Company  [Ref.  3], 
The  following  discussion  briefly  describes  the  basics  of  the 
flight  control  system  and  the  control  law  mechanization,  and 
is  intended  to  facilitate  understanding  the  model  development 

The  primary  flight  control  system  in  the  F/A-18  is  a 
fly-by-wire,  full  authority,  control  augmentation  system. 
The  control  law  computations  are  performed  by  four  flight 
control  computers  operating  in  parallel.   Each  computer  re- 
ceives input  from  the  aircraft  motion  sensors,  air  data 
computer,  and  pilot  stick  commands.   The  computer  operates 
on  the  input  signals  according  to  the  control  law  algorithms 
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and  outputs  the  command  signals  to  fly-by-wire  electron y- 
draulic  ser voac t ua t o r s .   Figure  2.1,  taken  from  Reference  3, 
shows  a  functional  block  diagram  of  the  flight  control  system 
Exclusive  of  angle  of  attack  and  air  data  sensors,  the  system 
has  quadraplex  redundancy.   The  system  provides  two  fail  op- 
erate performance  for  augmented  motion  feedback  control.   A 
third  failure  causes  the  system  to  revert  to  either  open  loop 
direct  electrical  link  control,  or  stabilator  mechanical 
control  . 

The  control  augmentation  system  is  gain  scheduled  with 
angle  of  attack  and  air  data  to  provide  optimum  flying  quali- 
ties throughout  the  flight  envelope.   Cross  axis  intercon- 
nects (e.g.,  rolling  surface  to  rudder  interconnect)  are 
provided  for  turn  coordination  and  maneuverability  at  high 
angles  of  attack.   The  control  system  also  provides  feedback 
to  counter  inertial  coupling  at  high  roll  rates. 

The  F/A-18  has  ten  primary  flight  control  surfaces: 
Right  and  left  stabilators,  leading  edge  flaps,  trailing 
edge  flaps,  ailerons,  and  rudders.   Longitudinal  control  is 
provided  by  collective  stabilator,  and  collective  leading 
and  trailing  edge  flaps.   Lateral-directional  control  is 
provided  by  differential  stabilator,  differential  leading 
and  trailing  edge  flaps,  ailerons,  and  rudders.   Collective 
leading  and  trailing  edge  flap  deflections  are  scheduled  by 
the  control  laws  and  are  a  function  of  angle  of  attack.   The 
flap  positions  are  designed  to  provide  optimum  L/D  during 
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Figure  2.1   Functional  Block  Diagram  of  Flight  Control  System 
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cruise,  and  improve  flying  characteristics  in  maneuvering 
and  high  angle  of  attack  flight.   Figure  2.2,  copied  from 
Reference  3,  shows  the  control  surface  positions  and  direc- 
tion for  positive  deflection. 

The  objective  of  the  thesis  did  not  require  a  complete 
model  of  the  F/A-18  control  system  as  given  in  Reference  3. 
The  following  assumptions  were  made  to  reduce  the  model 
complexity  . 

1)  The  aircraft  is  operating  in  the  up  and  away  flight 
phase.   Under  this  assumption  the  control  laws  are 
operating  in  the  a u t o- f la ps-up  configuration.   Con- 
trol law  configurations  for  the  takeoff  and  landing 
phases  were  not  modeled. 

2)  Only  inner  loop  control  is  modeled.   For  inner  loop 
control  the  pilot  provides  commands  to  the  system. 
The  auto  functions  (outer  loop  control)  were  not 
modeled  . 

3)  Control  is  provided  by  the  control  augmentation  system. 
The  unaugmented  modes  such  as  direct  electrical  link 
or  mechanical  backup  are  not  modeled. 

4)  The  failure  logic  provided  to  reconfigure  the  control 
laws  in  the  event  of  a  sensor  or  actuator  failure  is 
not  modeled. 

5)  The  aircraft  is  operating  with  gear  up,  speedbrakes 
in,  and  no  external  stores. 

6)  Spin  mode  control  logic  is  not  modeled. 

7)  The  aircraft  trim  system  is  not  modeled. 

8)  High  angle  of  attack  conditions  are  not  considered  in 
the  thesis  model,  therefore  control  law  configurations 
for  this  flight  condition  are  not  modeled.   For  this 
thesis  high  angle  of  attack  is  defined  as  flight  con- 
ditions above  15  degrees. 
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C.   F/A-18  DYNAMIC  MODEL  OVERVIEW 

The  functional  block  diagram  of  the  F/A-18  model  which 
couples  the  flight  control  system  to  the  basic  airframe  is 
shown  in  Figure  2.3.   The  diagram  represents  a  multi-input 
multi-output,  sampled  data,  closed  loop  control  system. 
Theory  on  the  analysis  of  sampled  data  systems,  as  the  one 
shown  in  Fig.  2.3,  is  extensive  and  covered  in  a  number  of 
texts  (see  references).   The  development  of  the  F / A -18  model 
assumes  the  reader  has  a  rudimentary  understanding  of  control 
theory,  and  in  particular  the  theory  of  sampled  data  control 
systems . 

In  the  nomenclature  used  to  represent  the  control  signals 
pilot  inputs  are  prefixed  by  the  letter  '  P  '  .   Actuating  sig- 
nals, and  signals  produced  by  the  aircraft  sensors  are  pre- 
fixed by  the  letter  'E'.   Control  surface  deflections  are 
prefixed  by  the  letter  'D'.   The  nomenclature  used  to  repre- 
sent the  aircraft  perturbed  motion  variables  and  control 
surfaces  is  given  in  Table  2.1.   To  denote  the  motion  axis 
which  is  being  controlled,  the  signal  will  be  suffixed  with 
either  x,  y,  or  z  to  denote  longitudinal,  lateral,  or  direc- 
tional axis  respectfully.   Finally  a  matrix  or  a  vector  will 
be  denoted  by  an  upper  case  letter.   A  scalar  will  be  denoted 
by  a  lower  case  letter. 

The  input  vector  shown  in  Fig.  2.3, 


P(t)t  =  [px(t)  py(t)  pz(t)] 


(2.1) 
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TABLE  2.1 
MOTION  VARIABLES  AND  CONTROL  SURFACE  NOMENCLATURE 


q  pitch  rate 

nz  normal  acceleration 

aa  angle  of  attack 

yr  yaw  rate 

rr  roll  rate 

ny  lateral  acceleration 

str  right  stabilator 

stl  left  stabilator 

ter  right  trailing  edge  flap 

tel  left  trailing  edge  flap 

ler  right  leading  edge  flap 

lei  left  leading  edge  flap 

ar  right  aileron 

al  left  aileron 

rr  right  rudder 

rl  left  rudder 


DEG/SEC 

G 

DEG 

DEG/SEC 

DEG/SEC 

G 

DEG 


represents  the  aircraft  longitudinal  and  lateral  stick,  and 
rudder  deflection  in  inches.   The  output  vector, 

Y(t)t=  [q(t)  nz(t)  aa(t)  yr(t)  rr(t)  ny(t)]   (2.2) 

represents  the  perturbed  motion  of  the  aircraft  about  some 
steady  state  operating  condition.   The  motion  variable  units 
are  degrees,  degrees/sec,  and  G's.   Each  block  in  the  dia- 
gram contains  a  mathematical  model  which  simulates  the  dynam- 
ics of  that  particular  component.   The  control  law  block 
contains  the  aircraft  flight  control  law  algorithms  modeled 
as  linear,  time  invariant,  discrete  state  equations.   One 
processing  channel  of  the  flight  control  computer  described 
above  is  represented  in  the  control  law  model.    The  input 
vectors  to  the  control  law  model  are  the  discrete  stick  and 
rudder  input  signals  in  inches, 
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P(k)t  =  [px(k)  py(k)  pz(k)] 


(2.3) 


nd  the  discrete  motion  feedback  signals  from  the  sensors, 


Es(k)t  =  [q(k)  nz(k)  aa(k)  yr(k)  rr(k)  ny(k)]    (2.4) 


(The  units  of  the  motion  feedback  signals  are  degrees,  degrees/ 
sec,  and  G's).   The  input  vector  to  the  control  law  equations 
is  therefore 


U(k)11  =  [Es(k)  |  P(k)] 


(2.5) 


The  output  vector  from  the  control  law  block, 

Ec(k)t  =  [estr(k)  estl(k)  eler(k)  elel(k)  eter(k) 

etel(k)  ear(k)  eal(k)  err(k)  erl(k)](2.6) 

represents  the  discrete  command  signals  to  the  flight  con- 
trol actuators  in  degrees.   The  actuator  command  signals  en- 
ter the  GAIN  block  which  represents  the  configuration  gain 
matrix  discussed  in  the  introduction.   The  state  space  equa- 
tions shown  in  the  actuator  block  model  the  dynamics  of  the 
flight  control  actuators.   The  input  vector  to  the  actuator 
block  , 

Ea(t)   -  [estr(t)  estl(t)  eler(t)  elel(t)  eter(t) 

etel(t)  ear(t)  eal(t)  err(t)  erl(t)]        (2.7) 

represents  the  continuous  time,  actuator  command  signals  in 
degrees.   The  output  vector  from  the  actuator  block 
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represents  the  control  surface  deflections  in  degrees, 

D(t)t  =  [dstr(t)  dstl(t)  dler(t)  dlel(t)  dter(t) 

dtel(t)  dar(t)  dal(t)  drr(t)  drl(t)]        (2.8) 

The  deflection  vector  is  input  to  the  airframe  small  pertur- 
bation model  represented  by  the  state  variable  equations  in 
the  airframe  block.   The  output  vector  from  the  small  per- 
turbation model,  Y(t),  enters  the  sensor  dynamics  block  which 
contains  the  state  variable  model  for  the  aircraft  rate  gyros, 
acceler ome t er s ,  and  angle  of  attack  sensors.   As  discussed 
above  the  sensors  output  the  feedback  signals  which  are  sent 
to  the  control  laws.   In  Chapter  III  the  mathematical  models 
which  simulate  the  dynamics  of  each  component  in  Fig.  2.3 
will  be  developed  in  detail. 

Analog  to  digital  converters  are  modeled  as  impulse  sam- 
plers.  It  is  assumed  that  all  samplers  are  operating  at  the 
same,  constant  sampling  rate.   (The  actual  system  uses  multi- 
rate  sampling.   In  the  thesis  model  only  a  single  sampling 
rate  is  used.   The  program  actually  allows  any  desired  sam- 
pling rate  to  be  input.)   The  mathematical  operation  per- 
formed by  the  impulse  samplers  is  shown  in  Figure  2.4. 
Modeling  the  analog  to  digital  converters  as  impulse  samplers 
is  a  valid  assumption  if  the  quantization  error  of  the  actual 
system  is  at  an  acceptable  level.   Digital  to  analog  conver- 
ters are  modeled  as  zero  order  hold  devices.   The  mathemati- 
cal operation  performed  by  the  zero  order  hold  is  shown  in 
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Fig.  2.5.   Finally  it  is  assumed  that  the  processing  time  be- 
tween the  sampler  inputs  and  zero  order  hold  output  is  very 
much  less  than  the  sampling  period  (i.e.,  no  processing  delay 
time  is  assumed).   This  assumption  is  used  in  transforming 
the  continuous  time  state  equations  (actuators,  airframe,  and 
sensors)  into  discrete  time  equations. 
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Figure  2.4   Mathematical  Operation  Performed  by  Impulse  Sampler 


-7'    0     /     21  37    4/    5/    o7     77     i 

Figure  2.5   Mathematical  Operation  Performed  by  Zero  Order  Hold 
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Ill 


MODEL  DEVELOPMENT 


A.   CONTROL  LAW  MODEL  DEVELOPMENT 

The  longitudinal,  lateral,  and  directional  control  law 
models  were  developed  from  the  block  diagrams  and  information 
provided  in  Reference  3.   In  addition  to  the  simplifying  as- 
sumptions listed  in  Chapter  II,  many  of  the  components  in 
this  system  were  eliminated  based  on  the  following 
considerations : 

1)  The  control  law  model  will  be  coupled  to  the  small 
perturbation  model  of  the  F/A-18.   It  is  assumed 
that  operation  of  the  control  system  will  remain 
within  a  linear  region.   Therefore  the  non-linear 
components  in  the  control  laws  were  eliminated. 
These  include  position  limiters,  rate  limiters, 
and  dead  band  regions.   Other  non-linear  functions 
in  the  system  which  were  essential  to  the  model 
(pitch  stick  gradient  for  example)  were  linearized 
by  a  taylor  series  expansion. 

2)  For  the  same  reasons  discussed  above,  the  portions 
of  the  system  which  provide  inertial  decoupling 
were  not  modeled. 

3)  In  the  model,  the  control  signals  are  input  as  dis- 
crete signals.  Therefore  stick  and  rudder  dynamics 
are  not  modeled. 

4)  Noise,  which  may  be  detrimental  to  control  system 
performance  due  to  aliasing,  is  not  introduced  into 
the  model.   Therefore  the  anti-aliasing  prefilters 
were  not  included  in  the  model. 

5)  The  structural  modes  were  not  included  in  the  F/A-18 
airframe  model.  Therefore  the  structural  notch  fil- 
ters were  not  included  in  the  model. 

6)  To  prevent  discontinuities  in  the  signals,  the  con- 
trol laws  utilize  faders  in  portions  of  the  system. 
Discontinuities  could  occur  during  start  up,  fail- 
ures, or  transitions.   Since  none  of  these  condi- 
tions are  included  in  the  thesis  model,  faders  have 
not  been  modeled. 
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7)   The  F/A-18  control  system  uses  multi-rate  sampling 
in  the  input  and  feedback  paths  (20,  40  and  80  hz 
sampling  rates  are  used).   To  develop  a  state  var- 
iable model  only  a  single  sampling  rate  was  con- 
sidered.  Therefore  the  iteration  averagers,  used 
to  mathematically  combine  two  discrete  signals  of 
different  sampling  periods  were  not  modeled.   The 
simulation  program  allows  any  desired  sampling  rate 
to  be  input.   For  this  thesis  80  hz  was  used  as  the 
sampling  rate. 

Figures  3.1  and  3.2  show  the  simplified  block  diagrams 
of  the  longitudinal  and  lateral-directional  control  laws. 
The  inputs  are  the  discrete  stick  and  rudder  commands,  and 
motion  feedback  signals.   The  output  signals  are  the  commands 
to  be  sent  to  the  flight  control  actuators  via  the  reconfig- 
urable  gain  matrix.   Together  Figs.  3.1  and  3.2  make  up  the 
control  law  block  shown  in  Fig.  2.3. 

The  blocks  in  the  control  law  model  represent  two  basic 
transfer  functions:   Function  gains,  and  digital  filters. 
Table  3.1  lists  the  notation  used  to  represent  the  transfer 
functions.   Function  oains  and  digital  filters  are  described 
below. 

. TABLE  3 .1 
CONTROL  LAW  TRANSFER  FUNCTION  NOTATION 
Prefix  Transfer  function 


Function  gain 

Longitudinal  digital  filter 
Lateral  digital  filter 
Directional  digital  filter 
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7.0 


-> 


F32a 
qc 


F12 
ri  ,  ps 


/v 


Longitudinal 
Input 


nz(k) 


-> 


P5 


Normal 
Acceleration 


F32a 
qc 

3.5 

i — '-r 

-^ 


F22 
qc 


F12 
ri  ,ps 


~*Z 


J-\ 


1 

q(k) 

P2 

Pitch    Rate 

F68 
qc 


F32a 

qc 


F12 
ri.ps 


P9 


-r 


F40 
qc , ps, ri 


->£ 


estx(k) 
— > 


A 


Collective 
Stabilator 


*& 


*% 


aa(k) 


^ 


Angle  of 
Attack 


Pll 


-> 


P12 


I.e.  flap  grad 


„ 1.328 


t.e.  flap  grad 


1.405 


elex (k) 


Collective 
Leading  Edge 
Flap 

etex (k) 

> 

Collective 
Trailing  Edge 
Flap 


Figure  3.1   Simpified  Longitudinal  Control  Law  Diagram 
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Figure  3.2   Simplified  Lateral-Directional  Control  Law  Diagram 
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1 .   Functions 

The  functions  in  the  control  law  diagrams  perform  the 
system  gain  scheduling  described  in  Sec.  II. B.  The  functions 
operate  on  the  air  data,  angle  of  attack,  and  normal  acceler- 
ation to  compute  the  system  gains.  Table  3.2  lists  the  nota- 
tion used  to  represent  the  function  inputs.  The  mathematical 
equations  which  define  the  gain  schedules  are  given  in  Refer- 
ence 3.  Copies  of  the  functions  used  in  the  control  law 
model  are  given  in  Appendix  A.  In  the  simulation  program  the 
function  gains  are  computed  using  steady  state  conditions  for 
all  input  values. 


ps 
qc 
r  i 
nz 
alpha 


TABLE  3.2 

FUNCTION  INPUT  VARIABLES 

Indicated  static  pressure  lbs/ft2 

Dynamic  pressure  lbs/ft2 

Pressure  ratio  (ps  /qc)  ND 

Normal  acceleration  G 

Angle  of  attack  DEG 


2 .   Digital  Filters 

Lead-lag  filters  are  used  in  the  system  to  shape  the 
output  response  and  provide  adequate  gain  and  phase  margins. 
An  integrator  is  used  in  the  forward  loop  of  the  longitudinal- 
system  to  provide  zero  steady  state  error  between  command  and 
feedback.   The  control  system  design  report  [Ref.  3]  gives 
the  filter's  continuous  time  transfer  function.   For  the  con- 
trol law  model,  the  digital  filter  coefficients  were  computed 
using  the  Tustin  transform.   All  filters  were  modeled  as  first 
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order  systems  with  two  numerator  coefficients  and  a  single 
denominator  coefficient.   For  example 


H(Z)  =  P9N1---Z  +  P9X2 
Z  -  P9D 


(3.1) 


is  the  model  for  filter  P9.   Appendix  B  gives  the  digital 
filters  used  in-  the  model,  and  the  method  used  to  compute  the 
coefficients . 

3  .   State  Space  Models 

Standard  control  system  analysis  techniques  were  used 
to  derive  the  state  space  models  for  longitudinal  and  lateral 
directional  control  laws  shown  in  Figs.  3.1  and  3.2.   The 
following  sections  outline  the  procedures  used. 

a.   Longitudinal  System 

The  block  diagram  in  Fig.  3.1  contains  five 
input/output  signal  flow  paths: 

1)  Pitch  rate  to  collective  stabilator 

2)  Normal  acceleration  to  collective  stabilator 

3)  Longitudinal  stick  to  collective  stabilator 

4)  Angle  of  attack  to  collective  leading  edge  flap 

5)  Angle  of  attack  to  collective  trailing  edge  flap 

To  obtain  the  individual  path  transfer  functions  the  signals 
are  mathematically  combined  to  give  the  following  three  Z  - 
transform  equations: 

Estx(Z)  =  Hl(Z)Eq(Z)  -  H2(Z)Enz(Z)  -  H 3 ( Z ) Epx ( Z )  (  3  .  2  ) 
Elex(Z)  =  H4(Z)Eaa(Z)  (3.3) 
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Etex(Z)  =  H5(Z)Eaa(Z) 


(3.4) 


The  expressions  for  each  transfer  function  are  given  in  Ap- 
pendix C.   To  obtain  a  state  space  expression  for  longitudi- 
nal control  laws,  the  individual  Z-transfer  functions  are 
first  expressed  in  state  variable  form.   The  state  variable 
equations  are  then  combined  according  to  equations  3.2  -  3.4 
This  procedure  is  as  follows: 

1)   The  transfer  function  for  the  pitch  rate  to  collective 
stabilator  path  can  be  expressed  as: 


Hl(z)  =  bOZ   +  blZ  +  b2 
(Z-P9D)  (Z-P2D) 


(3.5) 


The  numerator  coefficients  are  functions  of  the  system  gains 
and  filter  coefficients.   The  roots  in  the  denominator  are 
the  poles  from  filters  P9  and  P  2 .   (Appendix  C  gives  the  de- 
tailed expressions  for  the  numerator  coefficients). 

2)   The  state  space  representation  of  Eq.  3.5  is  obtained 
using  the  parallel  programming  method  outlined  in 
Reference  4. 


xl(k  +  l)" 
x2(k  +  2). 


P9D   0 
0   P2D 


'xl(k)- 

x2(k) 


eq(k)    (3.6) 


estxl(k)  =[qstl  q  s  t  2] 


'xl(k)- 
x2(k). 


+  qst3  eq ( k ) 


(3.7) 


Applying  the  same  procedures  to  the  remaining  Z-transfer 
functions  in  Eq .  3.2  results  in  similar  expressions: 
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x3(k+l) 
X4(k  +  2)J  = 


'P9D   0  " 
0   P5D 


x  3  (  k  )~ 

+ 

"1" 

x4(k)_ 

_1_ 

enz(k)     (3.8) 


estx2(k)  =  nzstl  nzst2 


'xl(k)" 
x2(k) 


+  nzst3  enz(k)    (3.9) 


for  the  normal  acceleration  path,  H2(Z);  and 

x-5(k  +  l)  =  P9D  x5(k)  +  1  px(k) 
estx3(k)  =  pxstl  x5(k)  +  pxst2  px(k) 


(3.10) 
(3.11) 


for  the  longitudinal  stick  path,  H3(Z). 

Appendix  C  details  the  procedures  used  to  com- 
pute the  coefficients  in  the  output  equations.   The  nomencla- 
ture used  to  represent  the  coefficients  in  the  output  equation 
combines  the  notation  of  the  input  signal  and  output 
control  surface,  followed  by  a  number  indicating  the  coeffi- 
cient's numerical  order  in  the  equation.   For  example: 

qstl,  qst2,  qst3 

are  the  first,  second,  and  third  coefficients  in  the  pitch 
rate  to  stabilator  output  equation  (Eq.  3.7).   With  this 
system  of  nomenclature  the  respective  signal  path  of  the  co- 
efficient is  easily  identified. 

3)   The  state  variable  equations  for  H1(Z),  H2(Z)  and  H3(Z) 
are  now  combined  according  to  Eq .  3.2. 

estx(k)  =  estxl(k)  -  estx2(k)  -  estx3(k)     (3.12) 


34 


Adding    Eqs.     3.7,     3.9    and    3.11    gives 


estx(k)    =  [qstl    qst2-nzstl-nzst2    -pxstl] 


x  1  (  k  )■ 

x2(k) 
x3(k) 
x4(k) 

x  5  (  k  ). 


Tqst3-nzst3] 


eq(k)  ' 
enz ( k ) 


+  -pxstl  px(k) 


(3.13) 


for  the  output  equation.   Equations  3.6,  3.8,  and  3.10  can 
be  combined  to  give 


xl(k+l) 

x2(k+l) 
x3(k+l) 
x4(k+l) 
x5(k  +  l)J 


P9D 


P2D 


P9D 


P5D 


P9D 


xl(k) 
x2(k) 
x3(k) 
x4(k) 
x5(k) 


"1 

0" 

1 

0 

0 

1 

0 

1 

_0 

0_ 

eq(k) 

enz ( k ) 


px(k) 


(3.14) 


The  state  space  representation  of  equation  3.2  is  given  by 
Eqs.  3.13  and  3.14.   These  equations  output  the  collective 
stabilator  command  given  the  longitudinal  stick  and  motion 
feedback  inputs.   Note  that  the  motion  feedback  and  longitu- 
dinal inputs  have  been  separated.   This  facilitates  coupling 
the  control  law  equations  to  the  aircraft  equations  to  be 
developed  later. 

4)   Similar  state  variable  equations  are  derived  for  the 
angle  of  attack  to  collective  flap  path  transfer 
functions  H4(Z),  and  H5(Z): 


x6(k  +  l)  =  PUD  x6(k)  +  1  eaa(k) 
elex(k)  =  aalel  x6(k)  +  aale2  eaa(k) 


(3.15) 
(3.16) 


for  the  AOA  to  collective  leading  edge  path  H4(Z);  and 
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x7(k+l)  =  P12D  x7(k)  +  1  eaa(k) 
etex(k)  =  aatel  x7(k)  +  aate2  eaa(k) 


(3.17) 
(3.18) 


for  the  A  0  A  to  collective  trailing  edge  path  H  5 ( Z )  . 

5)   The  state  equations  for  collective  stabilator,  and 
collective  leading  and  trailing  edge  flaps  are  now 
combined  to  give  the  following  state  variable  model 
for  the  longitudinal  control  laws: 


xl(k+l) 
x2(k+l) 
x3(k+l) 
x4(k+l) 
x5(k+l) 
x6(k+l) 
x7(k+l) 


P9D 


P2D 


P9D 


P5D 


P9D 


PUD 


P12D 


xl(k) 
x2(k) 
x3(k) 
x4(k) 
x5(k) 
x6(k) 
x7(k) 


1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

eq(k) 
enz ( k )  + 
eaa ( k ) 


px(k) 


(3.19) 


x  1  (  k  ) 

x2(k) 

estx  ( k  )~ 

"qstl 

qst  2 

-nzstl- 

nzs t  2 

-pxstl   0 

0  " 

x3(k) 

elex ( k  ) 

= 

0 

0 

0 

0 

0     aalel 

0 

x4(k) 

e  t  ex  (  k  )_ 

0 

0 

0 

0 

0      0 

aatel 

x5(k) 

x6(k) 

x7(k) 

"qst3-nzst3     0 
0    0     aale2 
0    0     a  a  t  e  2 


eq(k)  " 

enz ( k  ) 

+ 

eaa  (  k  )_ 

pxs  t  2 
0 
0 


px(k) 


(3.20) 


Equations  3.19  and  3.20  are  the  state  space  representation 
of  the  longitudinal  control  laws  shown  in  the  block  diagram 
in  Fig.  3.1.   The  equations  represent  the  longitudinal  con- 
trol law  model  which  computes  the  collective  stabilator 
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command,  and  collective  flap  commands,  given  the  longitudinal 
s tick  and  mo t ion  feedback  inputs. 

b.   Lateral-Directional  System 

The  procedures  outlined  above  are  applied  to  the 
Fig.  3.2  to  obtain  the  state  space  model  for  the  lateral- 
directional  system.  (To  complete  the  discussion  of  the  con- 
trol law  model  these  procedures  will  be  briefly  described.) 
The  individual  path  transfer  functions  are  first  obtained  by 
mathematically  combining  the  signals  in  Fig.  3.2  to  give  the 
following  equations: 

Esty(Z)  =  -H6(Z)Err(Z)  +  H7(Z)Py(Z)  +  H8(Z)Pz(Z)  (3.21) 
Eley(Z)  =  -H9(Z)Err(Z)  +  H10(Z)Py(z)       (3.22) 
Etey(Z)  =  -Hll(Z)Err(Z)  +  H12(Z)Py(z)      (3.23) 
Ea(z)  -  -H13(Z)Err(Z)  +  H14(Z)Py(Z)  +  H15(Z)Pz(Z)  (3.24) 
Er(Z)  =-H16(Z)Eyr(Z)  +  H17(Z)Err(Z)  + 
H18(Z)Eny(Z)  +  H19(Z)Py(Z)  +  H20(Z)Pz(Z)    (3.25) 

These  transfer  functions  represent  the  following  input/output 
signal,  paths,  numbered  respectfully: 

6)  Roll  rate  to  differential  stabilator 

7)  Lateral  stick  to  differential  stabilator 

8)  Rudder  pedal  to  differential  stabilator 

9)  Roll  rate  to  differential  leading  edge  flap 

10)  Lateral  stick  to  differential  leading  edge  flap 

11)  Roll  rate  to  differential  trailing  edge  flap 

12)  Lateral  stick  to  differential  trailing  edge  flap 
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13)  Roll  rate  to  aileron 

14)  Lateral  stick  to  aileron 

15)  Rudder  pedal  to  aileron 

16)  Yaw  rate  to  rudder 

17)  Roll  rate  to  rudder 

18)  Lateral  acceleration  to  rudder 

19)  Lateral  stick  to  rudder 

20)  Rudder  pedal  to  rudder 

The  expression  for  each  transfer  function  and  corresponding 
state  equation  are  given  in  Appendix  C.   Note  that  the  dif- 
ferential stabilator,  ailerons,  and  rudder  contain  the  trans- 
fer functions  for  the  cross  axis  interconnects  (e.g.,  H8(Z) 
and  H15(Z)  represent  the  rudder  to  rolling  surface  intercon- 
nect, and  H19(Z)  represents  the  rolling  surface  to  rudder 
interconnect.) 

The  state  equations  for  the  lateral  and  direction- 
al control  laws  are  given  as: 


"e  s  t  y  (  k  ) 

eley ( k  ) 

e  t  e  y  (  k  ) 

_ea(k)   J 


'0  -rrst  0' 

0  -rrle  0 

0  -rrte  0 

0  -rra  0 


eyr (k) 
er r ( k ) 
eny ( k ) 


'pyst  pzst 
py  le   0 
py  te   0 

pya  pza 


pv(k) 
pz(k) 


(3.26) 


for  the  lateral  system;  and 


x8(k  +  l)  ' 
x9(k+l) 
xl0(k+l) 
xll(k+l) 
.xl2(k  +  l). 


'Y3D 


Y3D 


Y5D 


Y5D 


Y5D 


r*8(k)  - 

x  9  (  k  ) 
x  1 0  (  k  ) 
xll(k) 
xl2(k) 
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[e  y  r  (  k  )" 
e  r  r  (  k  ) 
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"0 

0" 

0 

0 

0 

0 

+ 

1 

0 

_0 

1_ 

py(k) 

pz(k) 


(3.27) 


er(k)  =[-yrrl   rrrl   rrr2  pyrl  pzr^j 


x(8) 
x9(k) 
x  1 0  (  k  ) 
xll(k) 
x  1  2  (  k  ) 


[-  y  r  r  2 


r  r  r 


3  nyi] 


-ey  r  (  k  )■ 

e  r  r  (  k  ) 

_eny  (  k  ). 


py(k) 
+  [p  y  r  2  pzr2]        (3.28) 
pz(k) 


for  the  directional  system.   Note  that  the  lateral  system  is 
of  order  zero.   No  filters  were  included  in  the  model. 
c.   Combining  the  Control  Law  Models 

Equations  3.19  and  3.20,  and  3.26  -  3.28  are  com- 
bined to  give  the  3-axis  control  law  model. 


eq(k) 
enz ( k  ) 
12x1      12x12  12x1    12x6  eaa(k) 
Xc(k+1)  =    Ac   Xc(k)  +  Bfc  eyr(k)  + 

er r ( k ) 
eny ( k ) 


12x3  px(k) 
Be   py(k)   (3.29) 
pz(k) 


es tx ( k ) 
el  ex ( k  ) 
e  t  e  x  (  k  ) 
es  t  y ( k ) 
eley ( k ) 
e  tey ( k  ) 
ea(k) 
er(k) 


8x12  12x1     8x6 
Cc   Xc(k)  +  Dfc 


eq(k) 
enz ( k ) 
eaa ( k ) 
ey  r ( k ) 
er  r ( k ) 
eny ( k ) 


1x3  px(k) 

Be  py(k) 

pz(k) 


(3. 30) 


Equations  3.29  and  3.30  are  written  in  terms  of  the  matrix 
coefficients  in  Appendix  D.   These  equations  are  represented 
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by  the  discrete  state  equations  in  the  control  law  block  in 
Fig.  2.3. 

3  .   Variable  Gain  Matrix 

As  a  final  step  in  the  control  law  model  development, 
the  variable  gain  matrix  is  introduced.   The  command  signals 
which  are  output  from  the  control  laws  (Eqns.  3.29  and  3.30) 
are  distributed  to  the  right  and  left  actuators  according  to 
the  diagram  in  Fig.  3.3.   For  the  unimpaired  aircraft  the 
individual  gains  in  Fig.  3.3  are  set  to  unity.  ■  The  follow- 
ing matrix  gain  equation  represents  the  diagram  in  Fig.  3.3. 

'GAIN' 
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er(k) 


(3.31) 


The  matrix  gain  equation  will  be  recomputed  for  the  impaired 
aircraft  . 

B.   AIRCRAFT  MODEL  DEVELOPMENT 

A  3-axis  control  law  model  has  been  developed  which  will 
operate  on  the  stick  and  rudder  inputs  and  the  motion  feed- 
back signals.   The  model  outputs  the  control  signals  which 
are  distributed  to  the  aircraft  actuators  through  the  varia- 
ble gain  matrix.   The  next  step  in  building  the  F/A-18  model 
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Figure  3.3   Control  Law  Command  Signal  Distribution  and  Gains 
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is  to  develop  the  state  variable  equations  for  the  airframe, 
actuators,  and  sensors. 
1 .   Airframe  Model 

The  airframe  model,  obtained  from  the  Flight  Systems 
Branch  at  NATC,  is  a  linearized  small  perturbation  model  for 
both  the  longitudinal  and  lateral-directional  modes.   The 
dynamic  stability  and  control  derivatives  were  generated  from 
NATO's  F/A-18  Simulation  Package.   The  aircraft  state  for  the 
model  is  trimmed,  unacce ler a t ed  1-g  flight.   The  airframe 
state  equations  are  given  as: 


u(t) 
w(t) 
q(t) 

thed( t ) 

q(t) 
nz(t) 
aa(  t ) 


LONGITUDINAL 

4x4  u( t )  4x3  dstx( t ) 

Fx  w( t )  +   Gx . dlex( t ) 

q ( t )  dtex( t ) 
thed(t) 

3x4  u( t )  3x3  dstx(t) 

Hx  w(t)  +   Dx  dlex(t) 

q  (  t )  d  t  e  x  (  t  ) 
thed(t) 


(3.32) 


(3.33) 


LATERAL  DIRECTIONAL 
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v(t) 

4x5 

d  ley  ( 
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=  Fy  z 
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Gy  z 
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d  a  (  t ; 

phi(t) 

phi(t) 

dr(t; 
ds  t  y  ( 

t) 

r(t) 

3x4 

v(t) 

3x5 

dley  ( 

t) 

p(t) 

=  Hyz 

r(t) 

+ 

Dyz 

d  t  e  y  ( 

t) 

ny(t) 

p(t) 

da(t; 

phi(t) 

dr(t; 

(3.34) 


(3.35) 


Airframe  variable  definitions  and  units  are  listed  in  Table 
3.3.   The  definitions  of  the  stability  and  control  derivatives, 


Notation 


u(t) 

w(t) 
q(t) 
v(t) 
yr(t 
rr  ( t 
thed 
phi( 
dstx 
dlex 
d  tex 
ds  t  y 
d  1  e  y 
d  t  e  y 
da(t 
dr(t 
nz  (  t 
aa  ( t 
ny(t 


t) 
t) 
t) 
t) 
t) 
t) 
t) 
t) 


TABLE  3.3 
AIRFRAME  VARIABLE  DEFINITIONS 
Variable 


Units 


longitudinal  velocity  perturbation       ft/s 
perturbed  normal  velocity  ft/s 

perturbed  pitch  rate  rad/s 

perturbed  lateral  velocity  ft/s 

aircraft  perturbed  yaw  rate  rad/s 

aircraft  perturbed  roll  rate  rad/s 

aircraft  perturbed  pitch  angle  rad 

aircraft  perturbed  roll  angle 
collective   stabilator  deflection 
collective  leading  edge  flap  deflection 
collective  trailing  edge  flap  deflection 
differential  stabilator  deflection 
differential  leading  edge  deflection 
differential  trailing  edge  deflection 
aileron  deflection 
rudder  deflection 
normal  acceleration  ft/s2 

angle  of  attack  rad 

lateral  acceleration  ft/s2 


and  associated  units,  which  make  up  the  matricies  in  Eqs. 
3.32  through  3.35  are  given  in  Appendix  E.   Note  that  the 
units  of  the  perturbation  model  are  not  compatible  with  the 
actuators  or  sensors.   (The  actuator  output  units  are  in 
degrees,  and  the  sensor  input  units  are  in  degrees/sec,  de- 
grees, and  G's.)   The  input  and  output  variables  for  the 
perturbation  model  were  scaled  in  the  simulation  program  to 
properly  interface  the  models.   Combining  equations  3.32  - 
3.35  gives 
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u  (t) 

u(t) 

d  s  tx 

w(t) 

4x4 

4x4 

w(t) 

4x3  | 

4x5 

dlex 

q(t) 

Fx 

:   ° 

q(t) 

Gx 

0 

dt  ex 

t  h'e  d  (  t )   : 

thed(t) 

4- 

d  s  t  y 

T 

v(t) 

4x4 

4x4 

v(t) 

4x3 

4x5 

dley 

r(t) 

o 

Fyz 

y  r  (  t ) 

o      1 

Gy  z 

dtev 

p(t) 

p(t) 

da 

phi(t) 

• 

phi(t) 
u(t) 

dr 
ds  tx 

q(t) 

3x4  | 

3x4 

w(t) 

3x3 

3x5 

dlex 
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Hx 
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q(t) 

Dx 
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d  t ex 
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r(t) 

thed( t ) 
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j_ 

d  s  t  y 
dley 
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3x3 
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0 

Hyz 

r(t) 

0 

Dyz 

dtev 

ny(t) 

rr(t) 

phi(t) 

da 
dr 

(3.36) 


(3.37) 


To  study  r econ f igurab le  flight  controls  the  aircraft 
model  should  be  capable  of  using  the  full  set  of  control  sur- 
faces available  to  produce  the  required  forces  and  moments. 
To  achieve  this  the  control  surfaces  are  split  into  indepen- 
dent right  and  left  hand  complements  (i.e.,  right  elevator, 
left  elevator,  etc.).   The  equations  are  then  coupled  so  that 
a  complement  of  control  surfaces  used  either  collectively, 
differentially,  or  as  a  single  side,  will  produce  the  appro- 
priate moments.   For  example  the  stabilators  deflected  col- 
lectively will  produce  a  pitching  moment,  deflected 
differentially  will  produce  primarily  a  rolling  moment,  and  a 
single  side  deflected  will  produce,  to  some  degree,  moments 
about  all  three  axes. 

The  F/A-18  airframe  modeled  in  equations  3.36  and 
3.37  inherently  offers  control  coupling  through  the  stabila- 
tor,  leading  edge  flap,  and  trailing  edge  flap  surfaces. 
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Additional  longitudinal  coupling  could  be  achieved  with  the 
ailerons  and  rudders.   At  the  time  this  thesis  was  written 
control  derivatives  were  not  available  on  the  longitudinal 
effects  of  the  ailerons  or  rudders.   (The  rudder  is  capable 
of  toe-in  or  flare-out  and  will  effect  the  longitudinal  re- 
sponse of  the  aircraft.   This  feature  is  normally  used  during 
takeoff  and  landing.) 

To  split  the  control  surface  deflections  into  right 
and  left  hand  complements,  the  following  equations  are  used 
which  compute  the  deflection  inputs  to  Eqs.  3.36  and  3.37. 
(Also   refer  to  Fig.  2.2  which  shows  the  control  surface  posi 
tions  and  corresponding  positive  deflections.) 


LONGITUDINAL  DEFLECTIONS 

dstx  =.  (dstl  +  dstr  )  /  2 

dlex  =  (dlel  +  dler  )/2 

dtex  =  ( d tel  +  dter  )  /  2 


LATERAL-DIRECTIONAL  DEFLECTIONS 
dsty  =  dstl  -  dstr 
dley  =  -dlel  +  dler 
dtey  =  dtel  -  dter 
da  =  (drl  +  d  r  r  )  / 2 
dr  =  (dal  +  dar)/2 


Where  r  and  1  correspond  to  right  and  left  surfaces.   These 
equations  are  rewritten  in  the  following  matrix  format: 
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(3.38) 


'  LATD '  distribution  matrix 
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0 
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5_ 
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(3.39) 


Note  these  equations  are  for  the  unimpaired  aircraft  only! 
Damage  to  one  or  more  of  the  control  surfaces  will  change  the 
LONG  and  LATD  ma  trie ies  directly. 

Replacing  the  input  vectors  in  Eqs.  3.36  and  3.37 
with  the  r.h.s.  of  equations  3.38  and  3.39  gives  the  follow- 
ing modified  airframe  equations: 

Fm  Gm 

dstr 
u( t)  u( t  )  dstl 

w(t)       4x4    4 x 4  w ( t )       4x3   3x10  dler 
q(t)        Fx     0   q(t)        Gx   LONG  dlel  (3.40) 

th'ed(t)  = thed(t)  + dter 

v ( t )       4x4    4x4v(t)       4x5   5x10  dtel 
r(t)        0     Fyz  yr(t)      Gyz   LATD  dar 
p ( t )  rr(t)  dal 

phi ( t  )  phi ( t )  d  rr 

drl 
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Hm  Dm      d  s t  r 

u( t )  dstl 

q(t)       3x4  |  3x4  w(t)  3x3   3x10  dler 

nz(t)       Hx  |   0   q(t)  Dx   LONG  dlel 

aa(t)  =    thed(t)  + dter 

yr(t)      3x4    3x4  v(t)  3x5   5x10  dtel 

rr(t)       0     Hyz  yr(t)  Dyz   LATD  dar 

ny(t)                 rr(t)  dal 

phi( t  )  drr 

drl 


(3.41) 


The  names  of  the  individual  matricies  appear  above  the 
equations  . 

2 .   Actuator  Model 

Transfer  functions  for  the  flight  control  actuators 
were  given  in  Reference  3  and  are  listed  in  Appendix  F. 
These  transfer  functions  are  low  order  approximations  of  the 
more  complicated  actuator  models  presented  in  Reference  3. 
The  models  were  used  in  the  F/A-18  rigid  body  stability  anal- 
ysis described  in  Reference  3,  and  approximate  the  frequency 
characteristics  of  the  higher  order  models  out  to  5  hz. 

To  incorporate  the  actuators  into  the  modified  air- 
frame model,  the  transfer  functions  are  first  put  into  state 
variable  form.   This  procedure  is  outlined  in  Appendix  F. 
The  following  equations  represent  the  actuator  state  variable 
model : 


24x1    24x24  24x1    24x10  10x1 
Xa(t)  =  Fa  Xa(t)   +   Ga  Ea(t) 


10x1 
D(t) 


12x24 
Ha 


24x1 
Xa(t) 


(3.42) 


(3.43) 
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where , 

EaCt)11  =  [estr(t)  estl(t)  eler(t)  elel(t)  eter(t)  etel(t) 

ear(t)  eal(t)  err(t)  erl(t)] 

and  , 

D  (  t )  t  =     [dstr(t)  dstl(t)  dler(t)  dlel(t)  dter(t)  dtel(t) 

dar(t)  dal(t)  drr(t)  drl(t)] 

The  input  vector,  Ea(t),  represents  the  input  signal  in  de- 
grees from  the  control  laws,  via  the  GAIN  matrix.   The  output 
vector,  D(t),  represents  the  right  and  left  control  surface 
deflections  in  degrees.   Equations  3.40  -  3.43  are  combined 
to  give  the  following  equations  of  the  airframe  plus  actuator 
model  : 


8x1 
X  m  ( t  ) 

24x1 
X  a  (  t ) 


Fp 

8x8   |  8x10  10x24 
Fm      Gm    Ha 

24x8  |    24x24 
0         Fa 


Xp(t) 

8x1 
Xm(t) 

24x1 
Xa(t) 


Gp 

8x10 
0     10x1 

+   Ea(t) 

24x10 
Ga 


(3.44) 


q(t) 

nz(t) 

aa(  t 
yr(t 

rr(t) 
ny(t) 


Hp 


Xp(t) 


8x1 
6x8    6x12  12x24   Xm(t) 

Hm      Dm    Ha     

24x1 
Xa(t) 


(3.45) 
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As  before  the  names  of  the  matricies  appear  above  the  equa- 
tions.  The  airframe  plus  actuator  model  inputs  the  command 
signals  to  the  actuators,  and  outputs  aircraft  motion.   Since 
the  control  surface  deflections  in  the  perturbation  model  are 
in  radians,  it  was  necessary  to  scale  the  Hm  and  Dm  matricies 
to  interface  with  the  actuator  model  which  outputs  deflections 
in  degrees.   This  is  done  in  the  simulation  program. 
3 .   Sensor  Model 

In  the  final  step  of  the  development  of  the  aircraft 
model,  the  state  variable  model  for  the  aircraft  sensors  is 
incorporated  into  the  airframe  plus  actuator  model.   The 
sensor  transfer  functions  are  given  in  Appendix  G.   In  state 
variable  form  the  sensor  equations  are  given  as: 


q(t) 

nz(  t  ) 

llxl 

11x11 

llxl 

11x6 

aa(t) 

X's(t)   = 

Fs 

Xs(t) 

+ 

Gs 

r(t) 
p(t) 
ny(t) 

(3.46) 


eq(t) 
enz ( t ) 
eaa ( t ) 
eyr( t )  = 
er  r  ( t ) 
eny ( t  ) 


6x11   llxl 
Hs    Xs(t) 


(3.47) 


The  sensor  model  inputs  the  aircraft  motion  variables  in 
degrees/ second ,  degrees,  or  G's,  and  outputs  the  correspond 
ing  signals  to  the  control  laws  in  the  same  units.   Combin- 
ing Eqs.  3.44  -  3.47  gives  the  following  model  of  airframe 
plus  actuators  plus  sensors: 
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Fps 


Xps( t )     Gps 


32x1 

X'p  ( t ) 


llxl 
X  s  ( t ) 


32x32      32x11   32x1 
FP         0     Xp(t) 


11x6  6x32 
Gs    Hp 


11x11   llxl 
Fp    Xs(t) 


32x10 

Gp    10x1 

+   Ea(t) 

11x10 

0 


(3.48) 


eq(t) 
enz ( t ) 
eaa  ( t ) 
eyr  ( t ) 
er r ( t ) 
eny ( t ) 


Hps 

6x32  | 
0    ! 


Xps(t) 

32x1 
6x11  Xp(t) 

Hs    

llxl 
Xs(t) 


(3.49) 


The  aircraft  model  inputs  the  actuator  signals  from  the  con- 
trol laws  via  the  GAIN  matrix  and  outputs  the  motion  signals 
from  the  sensors  which  are  sent  to  the  control  law  .equations. 

C.   ASSEMBLING  THE  OVERALL  SYSTEM  MODEL 

A  mathematical  model  for  each  component  in  the  control 
system  shown  in  Fig.  2.3  has  now  been  developed.   Before  the 
individual  components  of  the  model  can  be  assembled,  the  dis- 
crete state  equations  for  aircraft  model  (Eqs.  3.48  and  3.49) 
must  be  computed.   Performing  this  operation  the  discrete 
state  equations  for  the  aircraft  are  given  as: 


43x1      43x43   43x1 
Xps(k+1)  =   Aps   Xps(k) 


43x10  10x1 
Bps   Ea(k) 


(3.50) 


6x1 
Es(k) 


6x44 
Hps 


43x1 
Xps(k) 


(3.51) 


Where 
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Xps(k) 
EaCk)1 


Es(k) 


[Xm(k) | Xa(k) | Xs(k) ] 
[estr(k)  estl(k)  eler(k)  elel(k)  eter(k)  etel(k) 

ear(k)  eal(k)  err(k)  erl(k)] 
[eq(K)  enz(k)  eaa(h)  eyr(k)  err(k)  eny(k)] 


The  Aps  and  Bps  discrete  matricies  are  computed  as  follows 


Aps  =  e 


F  p  s  ""  t  s 


F  p  s  "'"  s 
ps=    e      dsXGps 


(3.52) 
(3.53) 


Where  ts  represents  the  system  sampling  time. 

The  GAIN  matrix  is  now  introduced  to  interface  the  dis- 
crete aircraft  equations  with  the  control  law  equations.   Re- 
placing the  input  vector,  Ea(k),  in  the  aircraft  equation 
with  the  r.h.s.  of  the  GAIN  equation  (Eq.  3.31)  gives  the 
following : 


43x1      43x43   43x1 
Xps(k+1)  =   Aps   Xps(k) 


43x10  10x8   8x1 
+    Bps   GAIN  Ec(k) 


(3.52) 


6x1 
Es(k)  = 


6x44    43x1 
Hps   Xps(k) 


(3.53) 


For  convenience  the  control  law  equations  are  repeated 
below  : 


12x1      12x12  12x1       12x6   6x1       12x3  3x1 
Xc(k+1)  =    Ac   Xc(k)   +    Bfc  Es(k)   +    Be   P(k)   (3.29) 


8x1       8x12   12x1        8x6   6x1        8x3    3x1 
Ec(k)   =    Cc    Xc(k)   +    Dfc   Es(k)   +    Dc   P(k)  (3.30) 
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Where  : 

Xc(k-)t  =   [Xx(k)  |Xz(k)  ] 

Es(k)   =   [ e  q ( h )  enz(k)  eaa(h)  eyr(k)  err(k)  eny(k)] 
?(k)1    =       [px(k)  py(k)  pz(k)] 

Ec(k)t  =   [estx(k)  elex(k)  etex(k)  esty(k)  eley(k)  etey(k) 
ea(k)  er(k)] 

Equations  3.52,  3.53,  3.29,  and  3.30  are  now  combined  to 
give  the  following  matrix  equation: 

43x1       43x43  43x10  10x8  8x6  6x43  i  43x10  10x8  8x12  43x1 

Xps(k+1)     Aps  +  Bps   GAIN  Dfc  Hps   I    Bps  GAIN   Cc  Xps(k) 
_  i 

12x1  12x6  6x43         '      12x12        12x1 

Xc(k+1)  Bfc   Hps  I        Ac         Xc(k) 

43x10  10x8   8x3 
Bps   GAIN   Dc   px(k) 

+         py(k)  (3.54) 

12x3        pz(k) 
Be 

6x1        6x43     |    6x12  43x1        6x3 
Es(k)       Hps         0    Xps(k)       0   px(k) 
= +    ___  py(k)        (3.55) 

8x1      8x6  6x43     8x12   12x1        8x3  pz(k) 
Ec(k)     Dfc   Hps      Cc    Xc(k)        Dc 

These  equations  model  the  dynamic  response  of  the  F/A-18 
system  shown  in  Figure  2.3. 

D.   MODELING  EFFECTOR  IMPAIRMENT 

Effector  impairment  is  divided  into  four  groups  termed 
'effector  impairment  classes'  (EIC)  [Ref.  1].  Figure  3.4, 
copied  from  Reference  1,  defines  the  EIC's  and  indicates 
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STABILITY  & 
CONTROL  AFFECTED 


EIC  -   2 


E  IC  -  4 


Figure  3.4   Effector  Impairment  Classes 
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their  affect  on  aircraft  stability  and  control.   As  shown  in 
Fig.  3.4,  effector  impairment  class  one  is  unique  in  that 
only  aircraft  control  is  affected.   The  remaining  impairment 
classes  affect  aircraft  stability  as  well  as  control.   A 
self  repairing  control  system  must  be  capable  of  detecting 
and  classifying  effector  damage  to  compute  the  proper  recon- 
figuration gains. 

For  the  case  of  EIC=1,  one  or  more  of  the  aircraft  con- 
trol derivatives  will  be  altered.   In  the  F/A-18  model  de- 
veloped above  this  is  reflected  in  the  LONG  and  LATD  matricies 
(Eqs.  3.38  and  3.39).   For  example  if  the  right  stabilator  is 
impaired  the  elements  L0NG(1,1)  and  LATD(1,1)  would  be  set  to 
zer o  : 

LONG 


0 

.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  5 

.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  5 

.5 

0 

0 

0 

0 

LATD 


0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  5 

.5 

Compare  the  above  matricies  with  Eqs.  3.38  and  3.39  for  the 
undamaged  aircraft.   Note  that  the  damaged  system  will  now 
produce  a  lateral  input  for  a  given  longitudinal  command  and 
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vice  versa.   For  the  model  developed  in  this  thesis,  only 
class  one  effector  impairments  are  considered. 

E.   CONCLUSION 

In  this  chapter  the  mathematical  models  for  each  compo 
nent  in  the  block  diagram  of  Fig.  2.3  were  developed.  The 
individual  models  were  then  assembled  to  form  the  complete 
model  of  the  F/A-18  dynamical  system.  Next  the  simulation 
program  is  developed  to  compute  the  model  matricies  in  Eqs 
3.54  and  3.55,  and  compute  the  response  of  the  system  to 
stick  and  rudder  inputs. 
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IV.   PROGRAM  DEVELOPMENT  AND  MODEL  VALIDATION 

A.   INTRODUCTION 

To  validate  the  F/A-18  system  model  a  computer  program 
was  written  to  compose  the  model  matricies  in  Eqs.  3.54  and 
3.55,  and  compute  the  system  response  to  stick  and  rudder 
inputs.   In  addition  the  program  simulates  an  actuation  loss 
of  the  right  or  left  stabilator.   The  program  was  written  in 
VS  Fortran  on  the  IBM  3033  computer  at  the  Naval  Postgraduate 
School.   The  program  is  organized  to  offer  flexibility  for 
future  development  and  modification. 


B.   PROGRAM  STRUCTURE 

The  program  may  be  divided  into  four  major  operations: 

1)  Data  input 

2)  Air  data  computations 

3)  Composition  of  the  system  matricies 

4)  System  response  computation 

The  operations  are  performed  by  a  series  of  subroutines  as 
shown  in  the  flow  diagram  in  Fig.  4.1.   In  Fig.  4.1  the 
solid  lines  indicate  control  flow,  and  the  dotted  lines  indi- 
cate data  flow.   The  subroutines  perform  the  steps  outlined 
in  Chapter  II  to  compose  the  system  matricies  and  compute 
the  response.   A  functional  description  of  each  subroutine 
shown  in  Fig.  4.1  and  associated  variables  is  given  in  Appen- 
dix H.   Existing  subroutines  at  NPS  were  used  to  perform  the 
required  matrix  manipulations  [Ref.  5].   These  subroutines 
are  also  defined  in  Appendix  H. 
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FA18   DATA 

s. s. .flight 
conditions 

(stab.  &  cont. 
derivatives 


Interactive 

control  and 

failure 

parameters 
__ 

I 


FA13 

compose  the  total 
system  matricies 


7\ 


7K 


T 


FLITEl 

read  in  s . s . 
conditions  & 
cont.  &  failure 
parameters 


_*. 


I 


4-> 


LAWS 
compose  the 
control  law 
matrcies 


FLITE2 

read  in  a/c  stab 
and  cont.  derivs 


l 


GAIN 

compute  the 
variable  gain 
matrix 


A 


£ ^_ 


AIRDAT 

compute  air  date, 
cor  control  law 
input 


A 


± 


_  J 


EXPINT 

compute  discrete 
system  equations 


i 


-M- 


LONLAT 

compose  LONG  & 
LATD  matricies 

compose  unimpaired 
airframe  matrix 


k: I   I 


± 


ARCRFT 
compose  airframe  + 
actuator  +  sensor 
matricies 


A 


T 


MODEQ 

compose  modified 
airframe  matricies 


PLANT 

compose  airframe 
+  actuator 
matricies 


ACTU 

compose  actuator 
matricies 


I 


A 


SENSOR 

compose  sensor 
matricies 


- 


-  Data  Flow 
=  Control  Flow 


Figure  4.1   Simulation  Program  Flow  Diagram 
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Execution  of  the  program  is  controlled  by  the  'F18'  Exec 
program.   The  Exec  defines  the  input  and  output  data  files  as 
shown  in  Table  4.1.   No  user  options  are  provided  to  define 
data  input /output  files.   The  result  of  each  subroutine  com- 
putation is  written  to  the  FA18  RESULT  file,  and  can  be 
viewed  by  the  user  after  program  execution.   A  copy  of  the 
Exec  program  is  given  in  Appendix  I. 


Record 
Number 

01 


02 


03 


04 


TABLE  4.1 

INPUT/OUTPUT  FILE  DEFINITIONS 

Device    File       Type  Use 

Name 


Disk 


Disk 


Disk 


Disk 


FA18 


FA18 


DATA 


RESULT 


OPTMATD    DATA 


0PTPL0T    DATA 


Contains  the  F  A 1  < 
input  data 


program 


Contains  the  result  of 
the  operations  performed 
by  each  subroutine 

Contains  FA18  system 
matricies  (Eqs.  3.54  and 
3.55)  formated  for  the 
control  system  design 
package  at  NPS. 

Contains  the  system 
response  data  computed 
by  the  FA  18  program 
formated  for  the  inter- 
active plotting  program 
at  NPS. 


1  .   Data  Input 

Subroutines  'FLITE1'  and  'FLITE2'  perform  the  data 
input  operations.   All  data  is  read  from  the  FA18  data  file 
except  the  control  and  failure  parameters  which  are  input 
interactively  . 
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'FLITE1'  reads  in  the  aircraft  steady  state  flight 
conditions  and  control  and  failure  parameters.   As  indicated 
in  Fig.  4.1  the  flight  conditions  are  read  from  the  FA18 
DATA  file.   The  control  and  failure  parameters  are  read  in 
interactively.   This  allows  the  user  to  conveniently  run  the 
program  for  various  control  inputs  and  control  surface  fail- 
ures.   'FLITE2'  reads  in  the  basic  airframe  stability  and 
control  derivatives  from  the  FA  18  DATA  file.   The  stability 
and  control  derivatives  are  arranged  in  the  matrix  format 
shown  in  Eqs.  3.33  -  3.36,  and  in  Appendix  E. 

2  .   Air  Data  Computations 

Prior  to  composing  the  system  matricies,  air  data 
computations  are  performed  by  subroutine  'AIRDAT'.   The  sub- 
routine computes  the  air  data  inputs  to  the  control  law  func- 
tions using  the  standard  atmosphere  equations  [Ref.  6].   The 
program  does  not  include  the  logic  to  compute  atmospheric 
conditions  above  the  gradient  (troposphere)  region.   There- 
fore the  computations  are  valid  only  up  to  36,000  feet. 

3  .   Composing  the  System  Matricies 

As  seen  in  Fig.  4.1  the  system  matricies  are  com- 
posed by  ten  subroutine  operations.   The  operations  and 
associated  equations  in  Chapter  III  are  as  follows: 

1)  Compose  the  LONG  and  LATD  matricies  Eqs.  3.38  &    3.39 

2)  Compose  the  modified  airframe  matricies  Eqs.  3.40  & 
3.41 

3)  Compose  the  actuator  matricies  Eqs.  3.42  &    3.43 
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4)  Compose  the  sensor  matricies  Eqs.  3.46  &  3.47 

5)  Compose  the  airframe  plus  actuator  matricies  Eqs. 
3.44  &  3.45 

6)  Compose  the  airframe  plus  actuator  sensor  matricies 
Eqs.  3.48  &  3.49 

7)  Compute  the  discrete  system  matricies  Eqs.  3.50  &  3.51 

8)  Compose  the  GAIN  matrix  Eqs.  3.31 

9)  Compose  the  control  Law  matricies  Eqs.  3.29  &  3.30 
10)   Compose  the  total  system  maricies  Eqs.  3.54  &  3.55 

The  basic  component  matricies  are  composed  by  first  generat- 
ing a  null  matrix  of  proper  dimensions,  and  then  assigning 
the  coefficient  values  to  the  proper  elements.   These  matri- 
cies include: 

1)  LONG  and  LATD 

2 )  Actuator 

3 )  Sensor 

4)  Control  law 

5)  GAIN 

The  coefficient  values  to  actuator  and  sensor  models  are 
written  into  the  program  as  constants;  they  are  not  contained 
on  a  separate  data  file.   The  coefficients  in  the  LONG,  LATD, 
control  law,  and  GAIN  matricies  are  first  computed,  then 
assigned  to  the  appropriate  matrix  element. 

The  LONG  and  LATD  matricies  are  composed  by  subrou- 
tine 'LONLAT'.   As  explained  in  Section  III.C  these  matricies 
reflect  control  surface  damage  for  one  of  the  four  control 
surface  impairment  classes.   In  the  present  program  only 
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class  one  effector  impairments  are  simulated.   For  EIC=1  the 
appropriate  control  surface  coefficients  in  the  'LONG'  and 
'LATD'  matricies  are  set  to  zero.    'LONLAT'  contains  the  logic 
to  impair  the  right  or  left  stabilator.    'LONLAT'  also  com- 
poses the  unimpaired  airframe  control  matrix,  GmO .   This  ma- 
trix will  be  used  to  compute  the  impaired  gain  matrix  in 
subroutine  ' VGAIN'  . 

The  control  law  matricies  are  composed  by  subroutine 
'LAWS'.   The  subroutine  may  be  divided  into  four  operations: 

1)  Compute  function  gains 

2)  Compute  the  filter  coefficients 

3)  Compute  the  matrix  coefficients 

4)  Assign  coefficients  to  the  control  law  matricies  as 
in  Eqs.  D. 1  &  D. 2. 

The  gains  are  computed  according  to  the  function  definitions 

given  in  Appendix  A.   The  functions  use  the  air  data  computed 

in  'AIRDAT',  steady  state  AOA,  and  normal  acceleration  to 

compute  the  gains.   All  functions  are  programmed  exactly  as 

shown  in  Appendix  A  except  for  the  following  non-linear 

functions : 

Function  20  Longitudinal  stick  gradient 

Function  1  Lateral  stick  gradient 

Function  14  Rudder  pedal  gradient 

Function  42  RSRI  non-linear  gradient 

These  functions  were  programmed  using  the  linear  terms  of  a 
Taylor  Series  Expansion  about  the  origin.   Some  of  the  func- 
tions computed  in  the  program  are  not  used  in  the  control  law 
model.   They  are  provided  for  user  information  on  system  per- 
formance.  These  functions  are: 
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Function  24  Trailing  edge  flap  schedule 

Function  25  Trailing  edge  flap  schedule  qc  limit 

Function  27  Leading  edge  flap  schedule 

Function  29  Leading  edge  flap  schedule  ri  limit 

Function  37  Nz  limit  on  AOA  feedback 

Function  41  Rolling  surface  limit  schedule 

Function  112  Lateral  acceleration  gain 

Function  113  Lateral  acceleration  gain 

The  filter  coefficients  were  computed  using  the  procedures 

outlined  in  Appendix  B  to  transform  the  analog  filters  into 

digital  filters.   The  control  law  coefficients  were  computed 

using  the  procedures  outlined  in  Appendix  C  for  transforming 

the  control  path  transfer  functions  into  state  variable  form. 

The  function  and  coefficient  values  are  output  to  the  FA  18 

RESULT  file. 

Subroutine  '  V  G  A I N  '  composes  the-  GAIN  matrix  based  on 
control  effector  impairment.   For  the  unimpaired  system  the 
GAIN  matrix  appears  as  in  Eq .  3.31.   The  subroutine  is  de- 
signed to  implement  the  reconfiguration  algorithm  described 
in  Reference  1. 

Prior  to  composing  the  total  system  matricies  the 
discrete  form  of  the  continuous  system  matricies  must  be 
computed.   This  is  done  by  subroutine  'EXPINT'  [Ref.  5]. 
'EXPINT'  computes  the  matrix  exponential, 

F  p  s  *  t  s 


and  the  integral 


/   eFPs"sds 
o 
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The  discrete  system  equations  are  then  computed  as 


Aps  =  e 


F  p  s ";:"  t  s 


and 


t  s  i? 

u  s  F  p  s  "  s  ,     n 
,ps  =  f  e  ds  x  Gps 


The  Aps  and  Bps  matricies  are  written  to  the  OPTMATD  DATA 
file.   The  data  file  is  formatted  for  the  control  system 
design  package  at  NPS. 

4  .   Response  Computations 

System  response  is  computed  for  the  recursive  equa- 
tion 

55x1     55x55    55x1    55x3  3x1 
X(k+1)  =  AF18     X(k)  +  BF18  U(k) 

The  response  is  computed  for  500  data  points.   The  AF18  and 
BF18  matricies  are  defined  in  Eq.  3.54.   The  state  and  input 
vectors  are  defined  as: 


X(h) 


Xm(k)    Xa(k)    Xs(k)    Xc(k)  ] 


U(k)   =  [  px  py  pz  ] 

Response  data  for  all  55  states  is  written  into  the 
0PTPL0T  DATA  file  which  is  formatted  for  the  interactive 
plotting  routines  at  NPS.   The  variables  contained  in  the 
state  vector  which  are  relevant  to  the  user  for  viewing  are: 
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xl 

=  u 

x7  =  p 

x21 

= 

d  ter 

x2 

=  w 

x8  =  phi 

x23 

= 

dtel 

x3 

=  q 

x9  =  dstr 

x25 

= 

dar 

x4 

=  thed 

xl3  =  dstl 

x27 

= 

dal 

x5 

=   V 

xl7  =  dler 

x29 

= 

d  r  r 

x6 

=  r 

x!9  =  dlel 

x31 

= 

drl 

The  output  equation  (Eq.  3.55)  was  not  programmed. 

C.   PROGRAM  TESTING  AND  MODEL  VALIDATION 

To  test  the  simulation  program,  four  sets  of  runs  were 
made  for  the  following  control  inputs: 

1)  Positive  longitudinal  stick,  no  failure 

2)  Positive  longitudinal  stick,  right  stabilator  failed 

3)  Positive  lateral  stick,  no  failures 

4)  Positive  rudder  pedal,  no  failures 

All  deflections  were  -1  inch  step  inputs  of  3  second  duration 
All  runs  were  made  at  .6  mach/10000  feet.   The  sampling  rate 
was  set  to  80  hz.   (Additional  runs  were  made  at  40  hz  sam- 
pling rates  with  no  noticable  difference  in  the  output  re- 
sponse.)  Aircraft  steady  state  data,  including  the  stability 
and  control  derivatives,  are  given  in  Appendix  J.   The  FA18 
RESULTS  file,  including  the  function  and  control  matrix  co- 
efficient values,  is  given  in  Appendix  K. 
The  model  was  verified  for: 

1)  Correct  direction  of  motion  of  control  surfaces  and 
corresponding  aircraft  motion. 

2)  Expected  aircraft  response  for  a  right  stabilator 
failure . 

3)  Proper  augmented  aircraft  motion. 
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Reference  7  contains  response  plots  for  the  aircraft  for  simi- 
lar flight  conditions.   These  plots  were  used  to  verify  the 
model  for  the  proper  augmented  motion.   Some  of  the  plots 
were  reproduced  for  the  thesis  and  are  given  in  Appendix  L. 
The  control  inputs  used  in  Reference  7  are  slightly  different 
than  those  used  for  thesis  model.   Some  comparisons  however 
can  still  be  made.   All  response  plots  for  the  thesis  model 
are  given  in  Appendix  M. 

The  aircraft  response  to  a  positive  longitudinal  input 
for  a  no  failure  condition  is  shown  in  Figs.  M-l  -  M.3. 
The  motion  of  the  right  stabilator,  and  leading  and  trailing 
edge  flaps  compare  favorably  to  the  McDonnell  model  response 
shown  in  Fig.  L.l.   In  the  thesis  model  the  right  stabilator 
initially  travels  -2.5  degrees.   After  a  .5  second  transient 
period  the  stabilator  continues  to  -3.6  degrees  in  3  seconds. 
The  corresponding  McDonnell  response  shows  an  initial  travel 
of  -2.2  degrees.   After  the  .5  second  transient  period  the 
stabilator  continues  to  -1.8  degrees  in  3  seconds.   The 
major  difference  between  the  two  responses  is  the  stabilator 
rate,  the  thesis  model  being  slightly  greater  than  the 
McDonnell  model.   (Possible  explanations  for  the  discrepancies 
in  the  model  are  given  below.) 

The  leading  and  trailing  edge  flaps,  driven  by  angle  of 
attacks  feedback,  deflect  approximately  7.5  degrees  in  3.4 
seconds  as  shown  in  Figs.  M.l  and  M.2.   The  corresponding 
McDonnell  response  show  the  leading  and  trailing  edge  flaps 
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to  deflect  approximately  4  degrees  in  3.5  seconds.   (Note  the 
initial  conditions  of  the  flaps  show  Fig.  L.l  are  not  zero 
as  in  Figs.  M.l  and  M.2.)   Aircraft  pitch  rate  and  pitch 
angle  response  for  the  thesis  model  are  shown  in  Figs.  M.2 
and  M.3.   The  thesis  model  achieved  a  maximum  pitch  rate  of 
8  degrees/ second  compared  to  4  d eg rees / second  for  the 
McDonnell  model.   Also  the  McDonnell  model  achieved  a  con- 
stant pitch  rate  in  2.5  seconds.   The  thesis  model  achieved 
a  2.5  degree  pitch  attitude  in  5  seconds  compared  to  15 
degrees  shown  in. Fig.  1.1.   The  increased  travel  of  the 
stabilator  shown  in  Fig.  M.l  (compare  Fig.  L.l)  may  explain 
the  discrepancies  in  the  pitch  and  flap  responses. 

Figure  M.4  shows  the  response  of  the  aircraft  with  the 
right  stabilator  failed.   Note  the  decrease  in  pitch  rate 
and  pitch  angle  magnitudes.   Figure  M.4  shows  a  negative 
roll  rate  which  is  to  be  expected  for  a  right  stabilator 
failure.   Reference  7  did  not  give  aircraft  response  for 
control  surface  failures. 

Aircraft  response  to  a  positive  lateral  stick  input  is 
shown  in  Figs.  M.5  -  M.8.   The  motion  of  the  stabilator, 
aileron,  and  trailing  edge  flap  compare  favorably  to  the 
McDonnell  response  shown  in  Fig.  L.2.   The  aileron  in  the 
thesis  model  achieves  a  steady  state  deflection  of  1.1 
degrees  in  .5  seconds.   The  McDonnell  response  shows  a  steady 
state  deflection  of  1.6  degrees  in  .5  seconds.   The  thesis 
model  stabilator  (Fig.  M.5)  deflects  0.45  degrees  in  0.5 
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seconds,  the  corresponding  McDonnell  response  shows  a  0.6 
degree  deflection  in  0.5  seconds.   The  trailing  edge  flap 
in  Fig.  M.6  deflects  0.4  degrees  in  0.3  seconds  which  com- 
pares favorably  to  the  response  shown  in  Fig.  L.2.   (Leading 
edge  differential  flaps  are  not  used  at  the  flight  conditions 
tested.   Function  93,  which  sets  the  gain  on  the  leading  edge 
flap  path,  is  computed  as  zero.)   The  rudder  response  (Fig. 
M.6)  shows  an  initial  negative  deflection  which  is  expected 
for  coordinating  a  right  turn.   The  magnitude  of  the  deflec- 
tions for  the  two  models  is  approximately  the  same,  however 
the  thesis  response  is  slightly  more  oscillatory  than  the 
McDonnell  response.   The  aircraft  roll  rate,  yaw  rate,  and 
bank  angle  response  are  shown  in  Figs.  M.7  and  M.8.   Again 
only  slight  discrepancies  exist  between  the  two  models. 

The  aircraft  response  to  a  positive  rudder  pedal  input 
is  shown  in  Figs.  M.9  and  M.10.   The  corresponding  McDonnell 
response  is  shown  in  Fig.  L.3.   Comparison  of  the  two  models 
shows  the  shape  of  the  responses  to  be  approximately  the 
same.   However  the  magnitudes  in  the  thesis  model  are  very 
much  less  than  the  McDonnell  model.   For  example  the  rudder 
achieves  a  maximum  deflection  of  0.07  degrees  as  shown  in 
Fig.  M.9.   Compare  this  to  a  1.5  degree  deflection  for  the 
McDonnell  model  shown  in  Fig.  L.3.   The  response  of  the 
differential  aileron  and  stabilator  for  the  thesis  model 
(Figs.  M.9  and  M.10)  is  essentially  zero.   This  is  expected 
since  the  rudder  to  rolling  surface  interconnect  gain, 
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function  39,  is  computed  as  zero  for  low  angle  of  attack. 
(See  Appendix  A,  function  39.)   It  is  unknown  why  the  ailerons 
and  stabilators  respond  as  shown  in  Fig.  L.3.   Possibly  the 
response  shown  in  Fig.  L.3  is  based  on  a  different  flight 
control  program  than  that  used  in  the  thesis  model. 

Possible  explanations  for  the  discrepancies  in  the  re- 
sponse between  the  two  aircraft  models  are: 

1)  Improper  derivation  of  the  simplified  control  law 
model  based  on  the  assumptions  and  procedures  given 
in  Chapters  I  and  II. 

2)  Programming  discrepancies  resulting  in  erroneous 
computations . 

3)  Errors  in  the  aircraft  linear  small  perturbation 
model . 

A)   Differences  between  the  modeling  techniques  used 

in  Reference  7  and  in  this  thesis.   (The  MCAIR  model 
included  all  aerodynamic  and  control  system  non- 
linearities,  as  well  as  the  effects  due  to  digital 
time  delay  and  quantization.) 

5)  The  McDonnell  model  response  made  available  to  the 
author  may  be  based  on  a  slightly  different  flight 
program  than  used  in  the  thesis  model. 

The  author  investigated  each  of  the  items  listed  above.   It 
is  felt  the  assumptions  used  in  developing  the  control  law 
model  are  correct  given  the  available  information  on  the 
F/A-18  control  system.   The  simplifying  assumptions  used  in 
developing  the  model  may  account  for  some  of  the  differences. 
The  program  code  was  thoroughly  reviewed  and  revealed  no  dis- 
crepancies.  It  is  assumed  however  that  this  item  remains  a 
possible  source  of  error.   Two  possible  sources  of  error 
exist  with  respect  to  the  aircraft  perturbation  model: 
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1)  The  coefficients  in  the  NATC  small  perturbation 
model  are  computed  with  respect  to  the  aircraft 
body  axes  (see  Appendix  E).   This  may  effect  the 
magnitude  of  the  motion  variables  which  are  fed 
back  to  the  control  laws. 

2)  The  aircraft  model  developed  in  Reference  7  con- 
siders the  offset  position  of  the  ac eel er ome t er s 
from  the  aircraft  center  of  gravity.   It  has  not 
been  determined  if  this  effect  is  considered  in 
the  NATC  perturbation  model. 

Further  investigation  into  these  possibilities  should  be 

conducted  . 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

A  mathematical  model  was  developed  which  simulates  the 
dynamic  response  of  the  F/A-18  aircraft.   The  model  is 
designed  to  implement  the  reconfiguration  gain  matrix  for 
the  study  of  reconfigurable  control  systems.   A  program 
was  written  which  composes  the  discrete  time  state  variable 
matrices  of  the  aircraft,  and  computes  the  response  to 
stick  and  rudder  inputs.   The  program  also  simulates  the 
actuation  loss  of  the  right  and  left  stabilators.   Response 
plots  of  the  thesis  model  were  displayed  and  compared  to 
the  aircraft  model  developed  in  Reference  7.   Possible 
sources  of  error  were  discussed,  and  it  is  recommended  that 
further  investigation  into  each  of  these  areas  be  conducted. 

It  is  also  recommended  that  the  NATC  program  which  com- 
putes the  stability  and  control  derivatives  be  acquired  by 
the  Aeronautical  Engineering  Department  at  the  Naval  Post- 
graduate School,  to  be  made  an  integral  part  of  this  program 
This  would  expand  the  ability  of  the  program  to  simulate  the 
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response  of  the  aircraft  for  any  flight  condition,  and  various 
degrees  of  control  surface  damage. 
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APPENDIX    A 
FUNCTION    MATHEMATICAL    DESCRIPTIONS 
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FUNCTION    7 
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FUNCTION     10 
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FUNCTION     14 

RUDDER    PEDAL    GRADIENT 
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FUNCTION     17 
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FUNCTION    20 

LONGITUDINAL    STICK    GRADIENT 
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FUNCTION    24 
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FUNCTION    25 

TRAILING  EDGE  FLAP  SCHEDULE 
AUTO  FLAP  UP  (Qc  LIMIT) 


20 


LIMIT 
DEG 


10 


0 


0 


-<  (600.0,  17.0 


(G35.  0,5.0)  - 


500 


Or 


1000       1500 
•  LB/FT2 


2000 


C 


c 


PSF 


47.635  BIAS 


-0.05106 


84 


FUNCTION  27 
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.'   0.4 


RKF 


-i  -  I  .0  t 1 


■0.6913  BIAS 


4.0 


KFL 


LIMIT 
KFL 
0.05 


RKF  <  0 


—  K323 


RKF  >  0 
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FUNCTION    34 

DIFFERENTIAL    TRAILING    EDGE    FLAR 

GAIN    SCHEDULE 

AUTO    FLAP    UP     (AOA) 


Or  (O.O,      !   .0)  f- 


GAIN 
VOLTS/VOLT 


0.8 
0.6 
0.4 
0.2 


0.2 


-0.4^ 


2.5  BIAS 


AOA 
(DEG) 


-♦15,  1.0 


-*l  -5,  0.0) 


(10,  0.0)*- 


-4    -2 


AESOLUTE 
V*LUE 


2  4 

AOA  -DEG 


5  BIAS 


0. 


10 


GAIN 
VOLTS/VOLT' 
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FUNCTION    35 


LATERAL  FORWARD  LGGR  GAIN  SCHE 
AUTO  FLAP  UP  (AOA) 


)ULE 


GAIN 
VOLTS/VOLT 


)#0  10,  1.0! 


0.5 


114,  I .01 


7.7,  0.35) 


Z*>,    0.351 


tO        0         10        20 
AOA  -  DEG 


0.955  BIAS 


AOA 
IDEG1 


0.C325 


0.04221 


'• 


0.5  BIAS 


>  6  DEG 


r°" 


AOA  i    6  DEG 


LIMIT 

♦0.5 
♦0. 175 


30 


—I  2.0 


GA  IN 
I  VOLTS/VOLT  I 
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FUNCTION    35 

AILERON  GAIN  SCHEDULE 
AUTO  FLAP  UP  (Qc,  Ps ,  Rj ) 


Qr  (REVERSAL  POINT)  SCHEDULE 
1 1 .  a .  i ee  1 1 

~~ *l  2  .  S  ,     I  98  I  I 


2000r 


"y  10.7.     13331 

f^r,     'OCO-^-IO.O.     1333) 


HINGE  MOMENT  GAIN  SCHEDULE  IR,,  Pc 


co 


I  .0 


2.0  3.0 


10.43.  1 0 . 68 . 

1.01  I .01 

,of  <        » ps  i  4oo  pst 


GH 
I  VOLTS/VOLT  1 


0.5 


\        \ 
(0.8.    0.821 


10.8.    0.45! 


P<-    t    6*0    PST 


Pc    t    2116    PSF 


Ojr 


0.5 
R 


I  .0 


REVERSAL  GAIN  SCHEDULE  (Qc.  Gr) 

1.0  i < I  Or -200.     I.  01 


( Op ,    0,01 


Or -200  ucj 

0--     -     PSF 


^ [    SG-EC7     1  GA I N 

I  ulslMM    ;     ""  I  VOLTS/VOLT  1 


955    BIAS  Q^lPSF) 


L1W,T  I       ] 1       A       Or      A       F 

.0.7  ~ — i^  ^<y  i  o. 


IT 
10 

.0 


G» 


COOS 


7/27/BI        36/V6.0-8 
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FUNCTION    37 

NZ  LIMIT  ON  AOA  FEEDBACK 
AUTO  FLAP  UP  (NZA) 


1  ,5r 


GAIN 
VOLTS/VOLT) 


GAIN 
(VOLTS/VOLT) 
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FUNCTION    33 


RSP!    GAIN    SCHEDULE 
AUTO    FLAP    UP     IAOA.     Rj,    Pg ) 
O.ST 


l~».0.     0.2S7I 


-10 


•2.3X 


"1- 

I  LIMIT 

— |    ♦.5*9 

•.223 

rsxi 

rs  ■ 

aOaB 

tat 

• ,2«3oe  rs*i    ♦ 

•  -  10 

•  -.«3 

.33007 

AOA     ) 

aOam 

aoa« 

• .  i3B*«  rswi    - 

•  -S.4B? 

•  -.023 

.0177* 

AOA 
I0E3I 


LIMIT 

-3 


UPPER 

LIMIT 


)AOAM         i^5  .AO-ft 


*,  <0.3>» 


-0. 


»|>0.D4 
-1.2 


"  f^V-j" 


LIMIT 
U. 


*,  4  o.xyt 


AOA 


>    10   DEO  i 


h 


OAtN 

iva.rt^vaLTi 


0.O»«    BIAS  AOA    <     10    CEO 


I 


■0 


-f M 


0.72   BIAS 


— *MWh 


2.3    BIAS 
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FUNCTION  39 

RUDDER  PEDAL  TO  ROLL  CA5  INTERCONNECT 

GAIN  SCHEDULE 
AUTO  FLAP  UP  (AOA) 


.5 


.0 


GAIN 
VOLTS/VOLT 


AOA 

DEG 


0.5 


0 


0 


LIMIT 
+  25 
+  I  3 


/125, 


(13,  0.01/ 


0 


20 

AOA  -  DEG 


03329  EI  AS 


30 


40 


GA  I  N 

'0LT5/V0LT 
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FUNCTION    40 

PITCH    RATE    TEEDOACX    GAIN    SCHEDULE 

AUTO    FLAP    LP     ( Ps .     Oc ,    R  J  ) 


•AIM       0.*- 


Oc   j   «-»o  rsr   AMD  R|    <  0.73 


1SOO.     0.  1  I 


rs  -  rsr 


1.Z3    BIAS 


LIMIT 
0.73 


k? 


«i     <    0.73 


^ 


O.SS«=3    81  AS 


•  1 .3123    X    10- 


jt: 


LIMIT 

iaoo 

3O0 


LIMIT 

iaoo 
0.0 


-0.2BIT7    BIAS 


u 


•    B33    X     10"* 


AMD    «  |     1    0.  73  SOO    < 


r^ 


0.«73    8IAS 


— J    -«.S    A     10"* 


I .3Z3B    X     I0'3    BIAS 


-JOB    BIAS 


jfc-i-Rf 


•I  .  7l*2TI    X     10' 


>® 


S.r7*«    X     10"*    BIAS 


I  .  S»«l     X     10 


n— 6a 


F-^> 


0. 16SZ3    BIAS 


■3.H«I  3    X     |Q'S    I       l     < 


fk 


I .«T247    X     I0'J    BIAS 


-•.Z81  32    X     10' 


0^    <    «-*0 
OK   n,    <    0.73 


*io  i  soo 
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FUNCTION  4  1 

ROLLING  SURFACE  LIMIT  SCHEDULE 
AUTO  FLAP  UP  (AOA.Rj ) 


AOA_ 
IDEQf 


LIMIT 
IVOLTS) 


1.5 


1.0 


0.5 


(8,     1.01 


I  16.70,0. 10) 


122.    0.25)    R|{0.Z7 
119. 00. 0. 10)      R|^0.35 


0 

10 

20 
AOA  -  D£G 

30 

40 

B.O  BIAS 

1.0  BIAS 

LIMIT 
♦  22 

♦a 


0.052571 


L1UIT 

♦  0.48 

♦  0.27 


0.23533 


SV 


Hgh-S 


kJ 


ACA   ^       8.0 


-.00997 


0.5375 


LPFER 
LIMIT 
=  0.39 

1  .25 

.  / 

\ 

LOV£R 
LIMIT 


LIMI" 
(VOLTS 
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FUNCTION  42 

RSRI  NONLINEAR  GRADIENT 


8.  33.  3( 


Y  :  X  I  0. 0479 /X/+ 0.76 


-*-    X 


[- 18.33,  -30) 


0        20 
INPUT  -  VOLTS 


30 


-30-1- 


INPUT 
IVOLTS) 


LIMIT 
±18.33 


ABSOLUTE 
VALUE 


0.76  1  BIAS 


^^Zi_r-nKH[/ 


OUTPUT 
(VOLTS) 
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FUNCTION  45 

DIRECTIONAL  FORWARD  LOOP  GAIN  SCHEDULE 
AUTO  FLAP  UP  (Qq,  P5 ) 


1 .0 


GAIN 
I  VOLTS/VOLT 


0.5 


250       500 


Pc  1     1450  PSF 


I  1000,  0.5' 


750 


CCO  1250 


QC  -  PSF 


500  BIAS 


°C 
PSF) 


PS 
(PSF 


1.0  5  I  AS 


LIMIT 

*  l  450 

♦625 

0.00  12  121     | 

GA  !  N 
VOLTS/VOLT) 


hid 


FUNCTION    68 

PITCH  RATE  FEEDBACK  GAIN  SCHEDULE 
AUTO  FLAP  UP  (Qq) 


0.8r 

0.6h 
F68   0.4 
0.2 
0 


260,  0.655) 
,(274  KCAS) 


(480,  0.0) 
366  KCAS 


100  200  300  400  500 

Qr  -  PSF 


(PSF 


-480  BIAS 


F68 
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FUNCTION    90 

LATERAL  ACCELERATION  FEEDBACK 
GAIN  SCHEDULE  (Rj ,  Ps ) 


GAIN 
VOLTS/VOLT 


628    BIAS 

P*         .  ^    r- 

LIMIT 
♦366 
0.0 

fPSFl      "iO/ 

X 

n            -  l?V*\    m 

LIMIT 

♦2.08 

0.0 

1 

1 

1          vjy 

0.72    BIAS 

6.5    BIAS 


0.C2055I 


GAIN 
VOLTS/VOLT 
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FUNCTION    93 

DIFFERENTIAL  LEADING  EDGE  FLAP  GAIN  SCHEDULE 
AUTO  FLAP  UP  (Rj,  Psj.  NzAF 1 

(0.56.    0.06)  ^  ^10. 65,    0.061 


0.06 


0.04 


T33A 


0.02 


(0.434,    0.06 


(0.42,    01 
(0.294,    Ol/ 


10.74,    0.061 


0.2  0.4 


0.6 


O.B  I  .0 


I  .2 


•  Or 


F936      0.5  ■ 


IPSTI 


(7B5.     I .01  I I0O0,     I .0) 


1  .0 


(0.0,     1.0) 


F93C      0.5- 


'(628.    0.01 


5O0  I0O0 


Pe    -    PSF 


18.0,    0.51 


0.69    BIAS 


O.OOOIB 


429S7    — —I    *0. 

LJl 


LIMIT 

06 
0 


F93A 


628    6' AS 


LIMIT 
♦  157 
0.0 


6.3694X I0"3 


F93S 


A 


IFTfCW    FIGURE     I  I 


LIMIT 

♦5 

0.0 


0.  I 


F93C 


GAIN 
I  VOLTS/VOLT  1 


1.0  BIAS 
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FUNCTION     100 

SPEEDBRAKE    COMPENSATION     INCREMENT 
AUTO    FLAP    UP     (R j ) 


2.0 


I  .5 


DELTAHT 
(DEC)     *U 


0.5 


10.6,  2.0: 


10.4,  o.o; 


0.2       0.4 


4.0  BIAS 


0.6 
Ri 


-f  10.86,  2.0) 


10.96.  o.o: 


0.3       l  .0 


LIMIT 

0.96 

0.4 


m 

10.0 

-2C.0 

9.2  BIAS 


R|  i  0.7733 

^ 


Rj  >  0.7733 


LIMIT 
2.0 
0.0 


PI  00 
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FUNCTION  101 

differential  stabilizer  load  alleviation  schedule 
auto  Flap  up  (Rj,  ps,  n-af) 


ga:n 
i  volts/volt  i 


0.02 


0.01 


%F 


I  .0  G 


10.91  ,  0.0151, 
(0.79,  0.009), 


11.13. 
0.0151 


I  I .2.  0.015) 


0.2 


P5  t     1333  PST 


10.55.  O.Qi^ 
0.4       0.6 


\ 


Ps  >  176  1  PSF  ( ' -2.  0.OO9) 


^Ps  =  1000  PSF 
Ps     i     498  PSF 


o.a  i.o 


I  .2 


0.58   BIAS 

- 

(&\ 

LOWER 
LIMIT 
:    0.0 

n, 

{£} 

&H 

limit 

♦  935 
0.0 

— j\ 

71 

UPPER 
LIMIT 

:     >H01 

F>b       - 

n^ 

I 
498    BIAS 

798.5 

4  s:a 

s 

*-IOI 

A^ 

limit 

♦275.5 

♦  i65.3 

0.36949 

-&y- 

IFROU    FIGURF     1  I 

5. *4*6    X     lO"5 


HHh 


GAIN 

I  VOLTS/VOLT) 


IassoLjt 

VALUE 


l:m:t 

:     5.0 
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FUNCTION  96 

YAW  RATE  GAIN  SCHEDULE 
AUTO  FLAP  UP  IQc.  PS ' 


1.2 


I  134. 1.1)      PS  -    2<Z    1670,  1.1) 
I *-» -?+ 


100,     l-"\\        p5    ••    ^ 


1840. 1.1) 


GAIN 
(VOLTS/VOLT) 


0.8 


0.4 


uea.o.asi  isoo.o.asi 

1168.0.8),,      .     _„    ISOO.0.8) 


zoo     400     eoo     aoo     iooo    1200 
Qc  -  pst 

1 .54  BIAS 


UPPER 
LIMIT 
r     168 

A 

0.00*4  - 

r-r^ 

w 

o.ae  bia; 

1 

Uq    <    r>00    PSr 

UPPER 
LIMIT 

ipsri 

'  S — -C 

Qc    >    500    PSF 

1 

o.oooea 

Hg>- 

621 

3    BIAS 

LIMIT 
♦628 
♦242 

p^ 

—§9~ 

O.OOO 7772 

IPSM 

> * 

GAIN 
I  VOLTS/VOLT  I 
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FUNCTION  107 

LONGITUDINAL  I NERT I AL  GAIN  SCHEDULE 
AUTO  FLAP  UP  (Qc> 


0.0010 

1750.    0.0CO64) 

o.oooe 

/\\ 

GAIN 

0.00O6 

/                             \ 

(VOLTS/VOLT) 

O.CKXM< 

'10.0,    0.00041)     \ 

0.OO02 

\ 

°( 

(2220.    0.0 1\ 

)        400      8O0     1200    1600  20O0  2400 

Qc  -  pst 

750  BIAS 


8.4  X  10"4  SIAS 


IPST) 


ABSOLUTE 
VALUE 

-— 

5.714    X     10"7 

LOWER 
LIUIT 
-  0.0 


GAIN 
(VOLTS /VOLT) 
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FUNCTION     103 

DIRECTIONAL  I NERT I AL  GAIN  SCHEDULE 
AUTO  FLAP  UP  IQq) 


GAIN 
(VOLTS/VOLT) 


0.0l6r 


0.012 


0.008 


0.004 


(50,    0.0125) 


1390,    0.00185 


11200,       O.OOISS) 


J 


400  800 

QC    -    PSr 


1200 


0.01557    BIAS 


IPSF) 


LIMIT 
-390 

-50 

J 

X 

7.48    X     10*8 

-& 

1 

V- 

J 

i 

i 

D ♦ ^Jj     A      I U     " 

GAIN 
( VOLTS/VOLT  1 
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FUNCTION     114 

RUDDER  PEDAL  COMMAND  GAIN  INCREMENT 
AUTO  FLAP  UP  (Rj ) 


I5n 


GAIN 
(VOLTS/VOLT) 


.00 


LJUIT 
♦0.30 
♦0.  16 

.   f^h 

0.83286 

a.       ^MM 

"  toJ 

1  VOLTS/VOLT  1 

0.  1' 

1286 
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FUNCTION  1 13 


l.OO-, 


GAIN     .50H 
(VOLTS/VOLT) 


.00 


LATERAL  ACCELERATION  GAIN 
AUTO  FLAP  UP  IAOA) 


,oo 


— i — 

4.00 


i la.oo. i .00) 


a. 00  t2.00         18.00 

AOA-DEG 


1 1 

20. OO         24. 00 


AOA. 

IDEGI 


LIMIT 
+  18 
+  12 

(frl 

—  n  irrr7 

_   RA  IN 

) 

"1  U  .  I  OCjU  / 

y/sj 

IVOLTS/VCLT 

2 
BI 

0 
AS 
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APPENDIX  B 
DIGITAL  FILTER  MODELS 

A.   PITCH  RATE  LEAD-LAG  FILTER   P2 


.015(1+F22)S  +  1 

P2(S)  = — (P.l) 

.015S  +  1 


Using  the  tustin  transform  let 

2  (7-1) 

S  =   

ts(Z+l ) 

which  gives 


(P2N1)  (P2N2) 

2( .015) (1+F22)  +  ts        ts  -  2( .015) (1+F22) 
[ ]  Z  +  r ] 

2( .015)  +  ts                2(  .015)  +  ts 
P2(Z)  =   

2( .015)  -  ts 

z  _  r ] 

2( .015)  +  ts 


(B.2) 


(P2D) 


P.   NORMAL  ACCELERATION  LAG  FILTER   P5 


1 

P5(S)  =   (B-3) 

.04S  +  1 
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Using  the  Tustin  transform: 


(P5N1)       (P5N2) 

ts          ts 
1    + 

.OB  +  ts     .00  +  ts 
P5(Z)   =   (B.4) 

.00  -  ts 

z  -  r 1 

.00  +  ts 
(P5P) 


C.   INTEGRATOR   P9 


I 

P9(S)  =  

S 


(P. 5) 


Reference  3  approximated  this  filter  using  the  backward 
difference  method.   Let 


Z-l 

s  = 

tsZ 


which  gives 


(P9N1)  (P9N2) 

ts7  +  0.0 

P9(Z)  = (B.f) 

Z  -  1 

(PQP) 
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D.   LEADING  EDGE  AOA  LAG  FILTER   PI  1 


1 
Pll(S)  =   (B.8) 

0.39S  +  1 


Using  the  Tustin  transform 


(P11N1)  (P11N2) 

ts                ts 
f ]  z  + 

.39(2)  +  ts       .39(2)  +  ts 
Pll(Z)  =   (B.9) 

.39(2)  -  ts 

Z  -   

.39(2)  +  ts 

(PUP) 


E.   TRAILING  EDGE  FLAP  AOA  LAG  FILTER   PI  2 


1 
P12(S)  =  (B. 10) 

0.79S  +  1 


Using  the  tustin  transform 


(P12N1)  (P12N2) 

ts  ts 

r ]  z  + 

.79(2)  +  ts      .79(2)  +  ts 

P12(Z)  =   (B.ll) 

.70(2)  -  ts 
Z 

.79(2)  +  ts 
(P12N3) 
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F.   YAW  RATE  CANCELLER   Y3 


S 
Y3(S)  = (B.12) 

S  +  1 


Using  the  Tustin  transform: 


(Y3N1)  (Y3N2) 

2  -2 

r ]  z    +    

2  +  ts         2  +  ts 
Y3(Z)  =  

2  -  ts  (B. 13) 

Z   -   

2  +  ts 

(Y3P) 


G.     RSRI  LEAD  LAG  FILTER   Y5 


.75S  +  1 

Y5(S)  =   (B.14) 

.5S  +  1 


Using  the  Tustin  transform 


(Y5N1)  (Y5N2) 

2(  .75)  +  ts         ts  -  2( .75) 
[ "i  z   +   

2(  .5)  +  ts           2( .5)  +  ts 
Y5(Z)  =  

2( .5)  -  ts  (B. 15) 

Z 

2( .5)  +  ts 

(YSM3) 
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APPENDIX  C 
SIGNAL  PATH  TRANSFER  FUNCTIONS  AND  STATE  SPACE  MODELS 

A.   PITCH  RATE  TO  COLLECTIVE  STABILATOR  PATH  HI  (  Z  ) 

From  Fig.  3.1  the  forward  path  transfer  function  is 

H1(Z)    =    [P9*P2*F68*F32a*F12    +    P2* ( ( F68*F32a ) +F40 ) ]       (C.l) 

Let: 

A  =  F68*F32A*F12 

B  =  (F68*F32a)+F40 
Rewriting  Eg.  C.l  using  the  filter  Z-transform  expressions 
for  P9  and  P2  gives 


A*(P9N1*Z  +  P9N2)*(P2N1*Z  +  P2N2 ) 
H1(Z)  =   

(Z-P9D)*(Z-P2D) 


B*(P2N1*Z  +  P2N2) 

(C>2) 

Z-P2D 


Expanding  Eg.  C.2 

2 
H1(Z)  =  {[A*P9N1*P2N1  +  B*P2N1]  Z 

+   [A*(P9N1*P2N2+P9N2*P2N1 )  +  B* ( P2N2-P2N1 *P9D ) 1  Z 

+   [A*P9N2*P2N2  -  B*P2N2*P9Dl }  /  ( Z-P9D) * ( Z-P2D ) 

(C3) 
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Equation  C.3  can  be  expressed  in  terms  of  the  following 
partial  fraction  expansion: 


bO*Z  +  bl*Z  +  b2  qstl       ast2 

H1(Z)  = =  ast3  + + (C.4) 

(Z-P9D)*(Z-P2D)  (Z-P9D)    (Z-P2D) 


V/here 


qst3  =  bO 

2 
qstl  =  (bO*P9D  +  bl*(P9P)  +  b2 )  /  (P9P-P2P) 

2 
qst2  =  fbO*P2D  +  bl*fP2D)  +  b2 )  /  (P2D-P9P) 

Equation  C.4  can  now  be  put  in  the  following  state  variable 
form: 


[P9D   0  |xl(k) 
0   P2Djx2(k) 


(k) 


fC.5) 


estxl (k) 


=  g 


Stl  QSt 


+   qst3  q(k) 


(C.6) 


The  sytem  of  nomenclature  used  for  the  coefficients  in  the 
output  equation  ( Eq .  C.6)  is  described  in  Sec.    III. A. 

In  the  remainder  of  the  appendix  onlv  the  first  and  last 
steps  of  the  state  variable  development  will  be  given  for 
each  transfer  function.   The  equations  for  each  of  the 
coefficients  are  given  in  the  simulation  program  listing. 
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B.   NORMAL  ACCELERATION  TO  COLLECTIVE  STABILATOR  PATH  H2  (  Z  ) 

H2(Z)  =  [F12*F32a*3.5*P9*P2  +  3.5*F32a*P5]  (C.l) 


x3(k+l)      P9D   0   x3(k)      1 

+     nz(k)  (C8) 

x4(k+l )       0   P5D  x4(k)      1 


x3(k) 
estx2(k)   =  nzstl  nzst2         +   nzst3  cr  ( k )         (C.9) 

x4(k) 


PITCH  STICK  TO  COLLECTIVE  STABILATOR  PATH  H3 ( Z ) 

H3(Z)  =   r7.0*F32a*F12*P9  +  F20*F32a]  (CIO) 


Note:   F20  defines  the  longitudinal  stick  gradient.  The 
value  7.0  is  a  taylor  series  approximation 
about  the  origin  of  F20.  (See  F20  in  App.  A) 


x5(k+l)  =  P9D  x5(k)   +   1  px(k)  (C.ll) 

estx3(k)  =  pxstl  x5(k)   +   pxst2  px(k)  (C.12) 

D.    ANGLE  OF  ATTACK  TO  COLLECTIVE  LEADING  EDGE  FLAP  H4 ( Z ) 

H4(Z)  =  1  .328*P11  .  (C13) 


Note:  The  leading  edge  flap  schedule  is  defined  by  F27. 
The  value  used  in  Eg.  C.13  is  the  I.e.  flap 
gradient.   (See  F27  in  App.  A) 


x6(k  +  l)  =  PUD  x6(k)   +   1  aa(k)  (C.14) 

elex(k)  =  aalel  x6(k)   +   aale2  aa(k)  (C.15) 

117 


E.    ANGLE  OF  ATTACK  TO  COLLECTIVE  TRAILING  EDGE  FLAP  F5 ( Z ) 

H5(Z)  =  1.405*P12  (C.16) 


Note:   The  trailing  edge  flap  schedule  is  defined  by 

F24.  The  value  used  in  Eg.  C.16  is  the  t.e.  flap 
gradient.  (See  F24  in  App.  A) 


x7(k+l)  =  P12D  x7(k)   +   1  aa(k)  (C.17) 

etex(k)  =  aatel  x7(k)   +   aate2  aa(k)  (C.18) 

F.  ROLL  RATE  TO  DIFFERENTIAL  STABILATOR  PATH  H6 ( Z ) 

H6(Z)  =   rv7*(F4+rk6t)  (C.19) 

estyl(k)  =   rv7*(F4+rk6t)  rr(k)  (C.20) 

where : 

rv7  =  MIN(F6*F35,F6-F101 ) 

rk6t  =  LIMIT(0, . 12-F4, .12-F4) 
LL    UL      X 

G.  LATERAL  STICK  TO  DIFFERENTIAL  STABILATOR  PATH  H7 ( Z ) 

H7(Z)  =  3.22*rv7*F7*(F13+F4)  (C.21) 


Note:   The  lateral  stick  gradient  is  defined  by  Fl . 

The  value  3.22  is  a  taylor  series  approximation 
about  the  oriain  of  Fl .  (See  Fl  in  App.  A) 


esty2(k)   =  3 . 22*rv7*F7* ( F31+F4 )  py(k)  (C.22) 
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H.    RUDDER  PEDAL  TO  DIFFERENTIAL  STABILATOR  PATH  H8 ( Z ) 

H8(Z)  =  F14*1.33*rv7*F39  (C.23) 

esty3(k)  =  Fl 4*1 . 33*rv7*F39  pz(k)  (C.24) 

I.    ROLL  RATE  TO  DIFFERENTIAL  LEADING  EDGE  FLAP  PATH  H9 ( Z ) 

H9(Z)  =  F93*(F4+rk6T)  (C.25) 

eleyl(k)  =  F93*(F4+rk6t )  rr(k)  (C.261) 

J.    LATERAL  STICK  TO  DIFFERENTIAL  LEADING  EDGE  FLAP  H10(Z) 

H10(Z)  =  3.22*F93*F7*(F13+F4)  (C.27) 

eley2(k)  =  3 . 22*F93*F7* ( F31+F4 )  py(k)  (C.28) 

K.    ROLL  RATE  TO  DIFFERENTIAL  TRAILING  EDGE  FLAP  PATH  HI  1 ( Z  ) 
Hll(Z)  =  F31*F34*(F4+rk6t)  (C.29) 

eteylfk)  =  F31*F34* f F4+rk6t )  rrfk)  (C.30) 

L.    LATERAL  STICK  TO  DIFFERENTIAL  LEADING  EDGE  FLAP  HI  2 ( Z ) 

H12(Z)  =  3.22*F7*F31*F34*f F13+F4)  (C.31) 

etey2(k)  =  3 . 22*F7*F31*F34* f Fl 3+F4 )  py(k)  (C.321 
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M.    ROLL  RATE  TO  AILERON  PATH  HI  3 ( Z ) 

H13(Z)  =  F35*F36*.5*(F4+rk6t)  (C.33) 

eal(k)  =  F35*F36*.5*(F4  +  rk60  rrfk)  (C. 34) 

N.    LATERAL  STICK  TO  AILERON  PATH  H14(Z) 

H14(Z)  =  3.22*F7*F35*F36*.5*(F13+F4)  (C.35) 

ea2(k)  =  3.22*F7*F35  py(k)  (C.36) 

O.    RUDDER  PEDAL  TO  AILERON  PATH  H15(Z) 

H15(Z)  =  F14*F39*1 .33*F35*F36  fC.37) 

ea3(k)  =  F14*F39*1 . 33*F35*F36  pzfk)  (C.38) 

P.    YAW  RATE  TO  RUDDER  PATH  H16(Z) 

H16(Z)  =  F45*F96*COS(alpha)*Y3  (C.39) 

Note:   In  the  simulation  program  alpha  is  set  to  the 
steady  state  anqle  of  attack  in  degrees. 

X8(k+1)  =  Y3D  xP(k)  +  1  yr(k) 

erl(k)  =  yrrl  xR(k)  +  yrr2  yr(k) 
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0.    ROLL  RATE  TO  RUDDER  PATH  HI  7 ( Z ) 
H17(Z)  =  F45*F96*Y3*SIN(alpha) 

+  F45*F38*Y5*F30*0.76*f (2*rra)+(2*rrst) )   (C.40) 


Note:   The  second  term  on  the  r.h.s.  of  the  eauation 
results  from  the  rolling  surface  to  rudder 
interconnect  path.   The  value,  0.76,  represents 
F42,  the  RSRI  nonlinear  gradient.   It  is  a 
taylor  series  approximation  about  the  oriain 
of  F42. 


x2(k+l)     Y3D    0    x2(k)       1 

+       rrfk)  (C.41) 

x3(k+l)      0    Y5D   x3(k)       1 


x2(k) 
er2(k)  =   rrrl   rrr2  +   rrr3   rrfk)      (C.42^ 

x3(k) 


R.   LATERAL  ACCELERATION  TO  RUDDER  PATH  H18(Z) 

H18(Z)  =  F45*F90  (C.43) 

er3^k)  =  nyr  ny(k)  (C.44) 
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S.   LATERAL  STICK  TO  RUDDER  PATH   H19(Z) 

H19(Z)  =  F45*Y5*F38*F30*0.76*( (2*pya)+f 2*pyst) )   (C.45) 

X4(k+1)  =  Y5D  x4(k)   +   1  py(k)  (C.^6) 

er4(k)  =   pyrl  x4(k)  +  pyr2  py(k)  (C.47) 

T.   RUDDER  PEDAL  TO  RUDDER  PATH   H20(7) 

H20fZ)  =  F45*Y5*F38*0.76*F30*( (2*pza)+(2*pzst) 

-  F45*F14*( .5-(F17*F114) )       (C.48) 

Note:   The  first  term  on  the  r.h.s.  of  the  equation 
results  from  the  rudder  to  rollinq  surface 
interconnect  path. 

x5(k+l)  =   Y5D  x5(k)   +   1  pzfk) 
er5(k)  =  pzrl  x5(k)  +  pzr2  pzfk) 
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APPENDIX   D 


CONTROL  LAW  MATRICIES 


xl(k+l 

x2(k+l 

x3(k+l 

x  4  ( k + 1 

x5(k+l 

x6(k+l 

x7(k+l 

x8(k+l 

x9(k+l 

xlO(k+l) 

xll(k+l) 

xl2(k+l) 


[Ac] 


P9n 


P2D 


P9D 


P5D 


J 


P9D 


O 


PUD 


PI  2D1 


Y3D 


Y3D 


Y5D 


Y5D 


Y5P 


xl(k) 
x2(k) 
x3(k) 
x4(k) 

x  5  (  k  ) 

xo(k) 

x7(k) 

xR(k) 

x9(k) 

xlO(k) 

x  1  1  (  k  ) 

xl2(k) 


[Bfc] 


'Be' 


1 

0 

0  1 

I 

0 

0   ' 

0 

1 

o  1 

0 

1 

0 

1 

0 

q(k) 

0 

0 

o  1 

"" 

nz(k  ) 

0 

0 

1  1 

a  a  (  K  ) 

0 

0 

1 
-IT 

"o~ 

"  0 

yr(k) 
rr(k) 

I  o 

1 

0 

ny(k) 

0 

1  ° 

1 

0 

mt 

1  0 

0 

0 

I  o 

0 

0 

i! 

1 

0 

1 

n 

i 

1 

0 

1 

0 

1 

ro 

0 

I  ° 

** 

0 

i  1 

0 

0 
0 

0 
0       1 


px(k) 

PY  <  v-  ) 
nzfk) 


(d.i: 
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APPENDIX  E 
STABILITY  AND  CONTROL  DERIVATIVE  DEFINITIONS  AND  UNITS 


The  matricies  in  the  small  perturbation  model  (Eqs.  3.32 
through  3.36)  are  defined  in  this  appendix  along  with  the 
dimensions  of  the  stability  and  control  derivatives.   This 
information  was  obtained  from  the  flight  systems  branch  at 
the  Naval  Air  Test  Center. 


1 
XU 

zu 

1-ZWD 


FX 

MATRIX 

2 

3 

4 

XW 

XO-WB 

-G*CTHT 

zw 

(.-ZQ+UB) 
1-ZWD 

-G*STHT 

1-ZWD 

1-ZWD 

MWD*ZU        MWD*ZW        MWD* ( ZQ+UB )     MWD* ( -G ) *STHT 

3    MU+ MW+ M0+ 

1-ZWD  1-ZWD  1-ZWD  1-ZWD 


0 


0 


0 
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GX  MATRIX 


0  0  0 

ZDS  ZDLF  ZDTF 

1-ZWD  1-ZWP  1-ZWD 


MWD*ZDS  MWD*ZPLF  MWD*ZDTF 

3      MDS+ MDLF+ MDTF+ 

1-ZWD  1-ZWD  1-ZWD 


HX  MATRIX 


0 


zu 

zw 

zo 

1-ZWD 

1-ZWD 

1-ZWD 

0 

1 
UB 

0 
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DX  MATRIX 


0 


0 


0 


ZDS 
1-ZWD 


ZDLF 
1-ZWD 


ZDTF 
1-ZWD 


FYZ  MATRIX 


YV             YR-UB   .  YP+WR  G*CTHT 

1-YVD          1-YVD  1-YVD  1-YVD 

NVD*YV       NVD*( YR-UB)  NVD*(YP+WP)  NVD*G*CTHT 

2  NV+ NR+ NP  + 

1-YVD                           1-YVD  1-YVD  1-YVD 

LVD*YV       LVD* (YR-UB)  LVD*(YP+WB)  LVD*G*CTHT 

3  LV+ LR+ LP  + 

1-YVD          1-YVD  1-YVD  1-YVD 

STHT 

4  0              1  0 

CTHT 
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GYZ  MATRIX 


YDHT         YDLF        YDTF  YDA  YDR 

1-YVD         1-YVD       1-YVD  ]  -YVD  1-YVD 


NVD*YDHT     NVD*YDLF      NVD*YDTF     NVD*YDA      NVD*YDR 

NDHT+ NDLF+ NDTF+ NDA+ NDR+ 

1-YVD  1-YVD  1-YVD  1-YVD  1-YVD 


LVD*YDHT      LVD*YDLF      LVD*YDTF     LVD*YDA     LVD*YDR 

LDHT+ LDLF+ LDTF+ LDA+ LDP+ 

1-YVD         1-YVD  1-YVD  1-YVD         1-YVD 


HYZ  MATRIX 
12  3 

0  10 

0  0  1 

YV  YR  YP 

1-YVD         1-YVD      1-YVD 
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DYZ  MATRIX 

12        3        4  5 

0         0        0        0  0 

0         0        0        0  0 

YDHT     YDLF     YDTF     YDA  YDR 

1-YVD    1-YVD    1-YVD    1 - YVD  1-YVD 


The  following  abreviations  were  used  in  the  above 
matricies : 

UB  =  BODY  AXIS  LONGTTIDINAL  WIND 
WB  =  BODY  AXIS  VERTICAL  WIND 
CTHT  =  COS(THETA) 
STHT  =  SIN(THETA) 
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DFFINITION    OF  DIMENSIONAL  STABILITY  DERIVATIVES 


NAME  PARTIAL  DERIVATIVE  OF: 


xu 

LONGITUDINAL  FORCE 

xw 

LONGITUDINAL  FORCE 

XQ 

LONGITUDINAL  FORCE 

XWD 

LONGITUDINAL  FORCE 

XOSB 

LONGITUDINAL  FORCE 

XDTH 

LONGITUDINAL  FORCE 

XDS 

LONGITUDINAL  FORCE 

XDLF 

LONGITUDINAL  FORCE 

XDTF 

LONGITUDINAL  FORCE 

ZU 

VERTICAL  FORCE 

ZW 

VERTICAL  FORCE 

2Q 

VERTICAL  FORCE 

ZWD 

VERTICAL  FORCE 

ZDSB 

VERTICAL  FORCE 

ZDTH 

VERTICAL  FORCE 

ZDS 

VERTICAL  FORCE 

2DLF 

VERTICAL  FORCE 

ZDTF 

VERTICAL  FORCE 

nu 

PITCHING  MOMENT 

MW 

PITCHING  MOMENT 

MO. 

PITCHING  MOMENT 

MWD 

PITCHING  MOMENT 

MDS8 

PITCHING  MOMENT 

MDTH 

PITCHING  MOMENT 

MDS 

PITCHING  MCMENT 

MDLF 

PITCHING  MOMENT 

MDTF 

PITCHING  MOMENT 

YV 

LATERAL  FORCE 

YVD 

LATERAL  FORCE 

YR 

LATERAL  FORCE 

YP 

LATERAL  FCRCE 

YDA 

LATERAL  FORCE 

YDR 

LATERAL  FORCE 

YDLF 

LATERAL  FORCE 

YDHT 

LATERAL  FORCE 

YDTF 

LATERAL  FORCE 

LV 

ROLLING  MOMENT 

LVD 

ROLLING  MOMENT 

LR 

ROLLING  MOMENT 

LP 

ROLLING  MOMENT 

LDA 

ROLLING  MOMENT 

LDR 

ROLLING  MOMENT 

LDLF 

ROLLING  MOMENT 

LDHT 

ROLLING  MOMENT 

LDTF 

ROLLING  MOMENT 

WITH  RESPECT  TO; 

FORWARD  VELOCITY 
VERTICAL  VELOCITY 
PITCH  RATE 

VERTICAL  ACCELERATION 
SPEED  BRAKE 
THROTTLE 

HORIZONTAL  STABILATOR 
LEADING  EDGE  FLAPS 
TRAILING  EDGE  FLAPS 

FORWARD  VELOCITY 
VERTICAL  VELOCITY 
PITCH  RATE 

VERTICAL  ACCELERATION 
SPEED  BRAKE 
THROTTLE 

HORIZONTAL  STABILIZER 
LEADING  EDGE  FLAPS 
TRAILING  EDGE  FLAPS 

FORWARD  VELOCITY 
VERTICAL  VELOCITY 
PITCH  RATE 

VERTICAL  ACCELERATION 
SPEED  BRAKE 
THROTTLE 

HORIZONTAL  STABILIZER 
LEADING  EDGE  FLAPS 
TRAILING  EDGE  FLAPS 

LATERAL  VELOCITY 

LATERAL  ACCELERATION 

YAW  RATE 

ROLL  RATE 

AILERON 

RUDDER 

DIFFERENTIAL  LE  FLAPS 

DIFFERENTIAL  HORIZ  STABS 

DIFFERENTIAL  TE  FLAPS 

LATERAL  VELOCITY 

LATERAL  ACCELERATION 

YAW  RATE 

ROLL  RATE 

AILERON 

RUDDER 

DIFFERENTIAL  LE  FLAPS 

DIFFERENTIAL  HORIZ  STABS 

DIFFERENTIAL  TE  FLAPS 


14 


UNITS 

fT7SEC2 

FT/SEC2 

RAD/SEC2 

FT/SEC2 

RAD/SEC2 

PCT/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

FT/SEC2 

FT/SEC2 

RAD/SEC2 

FT/SEC2 

RAD/SEC2 

PCT/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

FT/SEC2 

FT/SEC2 

RAD/SEC2 

FT/SEC2 

RAD/SEC2 

PCT/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

FT/SEC2 

FT/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

FT/SEC2 

FT/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/SEC2 

RAD/5EC2 

RAD/SEC2 

RAD/SEC2 
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DEFINITION   OF  DIMENSIONAL  STABILITY  DERIVATIVES 

(Continued) 


NAME  PARTIAL  DERIVATIVE  OF: 


NV 

YAWING  MOMENT 

NVO 

YAWING  MOMENT 

NR 

YAWING  MOMENT 

NP 

YAWING  MOMENT 

NDA 

YAWING  MOMENT 

NDR 

YAWING  MOMENT 

NDLF 

YAWING  MOMENT 

NDHT 

YAWING  MOMENT 

NDTF 

YAWING  MOMENT 

WITH  RESPECT  TO: 


UNITS 


LATERAL  VELOCITY 

FT/SEC2 

LATERAL  ACCELERATION 

FT/SEC2 

YAW  RATE 

RAD/SEC2 

ROLL  RATE 

RAD/SEC2 

AILERON 

RAD/SEC2 

RUDDER 

RAD/SEC2 

DIFFERENTIAL  LE  FLAPS 

RAD/SEC2 

DIFFERENTIAL  HCRIZ  STABS 

RAD/SEC2 

DIFFERENTIAL  TE  FLAPS 

RAD/SEC2 
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APPENDIX  F 
ACTUATOR  TRANSFER  FUNCTIONS 


A.   STABILATOR 

S    2    2(0.068) 
( )   + s  +  1.0 

Dst(S)  82.0         R2.9 

Est(S)  S    2   2(0.41)  S   2    2(0.59) 

c( )   + s  +  li   r( )  + S  +  1? 

36.4       36.4  105.3      105.3 


The  procedures  outlined  in  Ref.  8  were  used  to  put  the 
.ransfer  functions  in  state  variable  form. 

1)   Arrange  the  transfer  function  as 

2 

b2  S  +  b3  S  +  b4 


4       3       2 
S  +  al  S  +  a2  S  +  a3  S  +  a4 


b2  =  2.1377E+03  a2  =  1.6122E+04 

b3  =  2.4101E+04  a3  =  4.9559E+05 

b4  =  1.4691E+07  a4  =  1.4691E+07 
al  =  1.5410E+02 

2)   Express  the  transfer  function  as  a  differential 
equation 

Dst(t)  +  al  Dst(t)  +  a2  Dst(t)  +  a3  Dst(t)  +  a4  Dst(t) 

b2  Est(t)  +  b3  Est(t)  +  b4  Est(t) 
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3)  Compute  the  state  equation  coefficients 

UO  =  bO  =  0 

Ul  =  bl  -  al*UO  =  0 

U2  =  b2  -  al*Ul  -  a2*U0  =  2.1377E+03 

U3  =  b3  -  al*U2  -  a2*Ul  -  a3*U0  =  -3.0532E+05 

U4  =  b4  -  al*U3  -  a2*B2  -  a3*Bl  -  a4*B0  =  2.7277E+07 

4)  Arrange  in  the  following  state  variable  form 

Fst  Gst 

xl  0    1  0  0  xl      Ul 

x2  =    0  0  1  0  x2  +   U2   Estft) 

x3  0  0  0  1  x3      U3 

x4  -a4  -a3  -a2  -al  x4      U4 

Hst 

xl 
Dst(t)  =10  0  0   x2 

x3 
x4 


Note  input  and  output  for  all  actuator  models  is  in 
degrees.   The  remaining  state  variable  models  are  developed 
as  above. 
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LEADING  EDGE  FLAP 


Dle(S)  1.0 


Ele(S)  S  S 

( +  1.0)  ( +  1.0) 

26.9  82.9 


xl(t)         xl(t) 

=   Fie        +  Gle  Ele(t) 
x2(t)         x2Ct) 


xl(t) 

Dle(t)  =  Hie 

x2(t) 


TRAILING  EDGE  FLAP 


DtefS)  1.0 


Ete(S)        S    2     2(0.71) 

( )   +   S   +   1  .0 

35.0        35.0 


xl(t)        xl(t) 

=  Fte        +  Gte  Ete(t) 
x2(t)         x2(t) 


xl(t) 

Dte(t)  =  Hte 

x2(t) 
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D.   AILERON 


Da ( S )  1.0 

Ea(S)        S    2  2(0.59) 

( )   +   S  +  1  .0 

75.0  75.0 


xl(t)       xlft) 

=  Fa        +  Ga  Ea ( t ) 
x2(t)       x2(t) 


xlft) 
Da(t)  =  Ha 

x2(t) 


E.   RUDDER 


Dr(S)  1.0 

Er(S)  S   2  2(0.69) 

( )   +  S  +  1.0 

72.1  72.1 


xl(t)       xl(t) 

=  Fr        +  Gr  Er(t) 
x2(t)       x2(t) 


xl(t) 
Dr(t)  =  Hr 

x2(t) 
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The  actuator  state  variable  model  natricies  in  ^os.  3.43 
&  3.44  are  arranged  as  follows: 

FA  MATRIX 


Fst 

Fst 

Fie  0 

Fie 

Fte 

Fte 
0  Fa 

Fa 

Fr 

Fr 


GA  MATRIX 


Gst 

Gst 

Gle  0 

Gle 

Gte 

Gte 
Ga 
0  Ga 

Gr 

Gr 


HA  MATRIX 


Hst 

Hst 

Hie 

Hte  0 

Hte 

Hie 

Hie 
0  Ha 

Ha 

Hr 

Hr 
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APPENDIX  G 
AIRCRAFT  SENSOR  TRANSFER  FUNCTIONS 


A.   RATE  GYRO  TRANSFER  FUNCTION 


S 

+  1.0 

Y(S)  131.7 

E(S)  S  S 

[ +  i.o]  r +  l.oi 

167.8  461.7      • 


B.   ACCELEROMETER  TRANSFER  FUNCTION 


S 

+    1.0 

Y(s)  235.8 

E(S)  S         2  2(.96) 

( )       +      S    +    1.0 

395.3  395.3 


ANGLE  OF  ATTACK  SENSOR  TRANSFER  FUNCTION 


Y(S)         1.0 

E(S)      S 

+  1.0 

14.0 


Note:  Transfer  functions  input  degrees/second  for  degrees) 
and  output  degrees/second  (or  degrees).   The  state  variable 
model  for  the  transfer  functions  in  Eqs.  3.46  &  3.47  are 
developed  using  the  procedures  outlined  in  appendix  F. 
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APPENDIX  H 
SIMULATION  PROGRAM  SUBROUTINES 


Subroutine:   FLITE1 


Description:   Reads  in  steady  state  fliqht  conditions,  stick 
and  rudder  commands,  and  control  surface 
failure  parameters. 


Calling  sequence:   CALL  FLITE1 (mach, alt , alpha , nz , ncont , nst , 

nstp, amp, rstf, lstf ) 


Input  arguements:  None 

Output  arguements: 

mach  Steady  state  mach  number 

alt  Steady  state  altitude  in  feet 

alpha         Steady  state  AOA  in  degrees 

nz  Steady  state  normal  acceleration.   Normally  1.0 

ts  Sampling  time  in  seconds 

ncont         Control  number:    l=Longi tud inal  stick 

2=Lateral  stick 
3=Pudder  deflection 

nst,  nstp     Control  start  time,  and  control  stop  time. 

Nst  and  nstp  are  actually  the  iteration  number: 
Start  time  fsecs. )=nst*ts 
Stop  time  ( sees . )=nstp*ts 

amp  Control  amplitude  given  as  deflection  in  inches 


rstf,  lstf    Right  and  left  stabilator  failure  parameters 
respectfully:   No  failure  =  0 

Failure  =  1 
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Subroutine:   FLITE2 


Operation:   Reads  in  basic  airframe  matricies  from  the 
'F/A-18'  data  file. 

Calling  arguments:   CALL  FLITE2 ( Fx , Gx , Hx , Dx , Fyz , Gyz , Hyz , Dyz , 

Nf x , Ngx , Nhx , Ndx , Nf yz ; Ngyz , Nhyz , Ndyz ) 


Input  arguments:   None 

Output  arguments: 

Fx ,  Gx,  Hx ,  Dx  Longitudinal  state  matricies 

Fyz,  Gyz,  Hyz,  Dyz        Lateral-directional  state  matricies 


Nf x, Ngx, Nhx , Ndx  Two-dimensional  vectors  giving  the 

Nfyz, Ngyz , Nhyz , Ndyz       number  of  rows  and  columns  of  the 

respective  matricies.  Example, 
Nfx(l)  =  Number  rows  in  Fx 
Nf x ( 2 )  =  Number  columns  in  Fx 
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Subroutine:   LONLAT 


Description:   Generates  the  LONG  and  LATD  matricies 

(Eqs.  3.38  and  3.39)  based  on  the  control 
surface  failure  parameters. 


Computes  the  unimpaired  input  matrix,  GmO,  to 
the  modified  ai frame  equations.   The  GmO  matrix 
will  be  used  in  the  reconfiguration  algorithm. 


Calling  Arguments:    CALL  LONLAT( r stf , Istf , i f ai 1 , Gx , Gvz , LONG, 

LATD, GmO , Nlong, Nlatd , NgmO , Ngx , Ngyz ) 


Input  arguments: 

rstf,  lstf       Right  and  left  stabilator  failure  parameters 

respectfully.   (See  subroutine  FLITE1 ) 

ifail  Failure  flag:   No  failure  =  0 

Failure  =  1 

Gx,  Gyz  Longitudinal  and  lateral-directional  input 

matricies . 

Ngx,  Ngyz        Two  dimensional  vectors  giving  the  number  of 

rows  and  columns  of  the  Gx,  and,  Gyz 
matricies . 


Output  arguments: 

LONG,  LATD       Matricies  described  in  Section  III.B.  which 

split  the  control  surface  deflections  into 
right  and  left  hand  sides.   The  LONG  and  LATD 
matricies  reflect  control  surface  failure  or 
damaqe. 

GmO  The  unimpaired  input  matrix  to  the  modified 

airframe  equations  ( Eq .  3.40).   GmO  is 
composed  usinq  the  unimpaired  LONG  and  LATD 
matricies.   It  is  subsequently  used  in  the 
reconfiguration  alaorithm. 

Nlong,  Nlatd     Two  dimensional  vectors  givino  the  number  of 
NgmO  rows  and  columns  in  the  respective  matricies. 
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Subroutine:   AIRDAT 


Description:   Computes  the  air  data  inputs  to  the  control 
laws  based  on  standard  atmosphere  conditions 


Calling  arguments:   CALL  AIRDATfmach, alt , temp, rho, ps i , a , 

ri) 


Input  arguements: 

mach  Steady  state  mach  number 

alt  Steady  state  altitude  in  feet 

Output  arguements: 

temp  Standard  atmosphere  temperature  in  degrees 

rankine 

rho  Standard  atmosphere  densitv  in 

lb-sec**2  /  ft**4 

psi  Standard  atmosphere  pressure  in  lb/ft**2 

a  Sonic  velocity  in  ft/sec 

qc  Dynamic  pressure  in  lb/ft**2 

ri  Pressure  ratio  =  gc/psi 
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Subroutine:   MODEQ 


Description:   Composes  the  modifier?  airframe  matricies  (Eos 
3.40  and  3.41) 


Calling  Arguments:   CALL  MODEQ ( Fx , Gx , Hx , Dx, Fyz , Gyz , Hyz , Dyz , 

LONG,  LATD,  Fm,Gm,  Hm,  Dm,  Nf  x  ,  Ngx  ,  Nhx  ,  Ndx  , 
Nf yz , Ngyz , Nhyz , Ndyz ) 


Input  arguments: 

Fx,  Gx,  Hx,  Dx 
Fyz,  Gyz,  Hyz,  Dyz 


LONG ,  LATD 


Longidudinal  and  lateral-directional 
basic  airframe  matrici.es  (Fas.  3.32 
-  3.35)   Read  in  from  subroutine 
FLIGHT1 

LONG  and  LATD  matricies  (Eas.  3.38 
and  3.39)  Generated  in  subroutine 
LONLAT 


Nfx,  Ngx,  Nhx,  Ndx        Two  dimensional  vectors  giving  the 
Nfyz,  Ngyz,  Nhyz,  Ndyz    number  of  rows  and  columns  in  the 
Nlong,  Nlatd,  .respective  matrix. 


Output  arguements: 
Fm,  Gm,  Hm,  Dm 

Nfm,  Ngm,   Nhm,  Ndm 


Modified  basic  airframe  matricies 
(Eqs.  3.40  and  3 .41 ) 

Row  and  column  vectors  as  described 
as  described  above 
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Subroutine:   SENSOR 


Description:   Composes  the  sensor  matricies  (Eqs.  3.47  and 
3.46) 


Calling  arguments:   CALL  SENSOR( Fs , Gs , Hs , Nf s , Ngs , Nhs ) 

Input  arguments:   None 

Output  arguements: 

Fs ,  Gs,  Hs        Sensor  matricies 

Nfs,  Ngs,  Nhs      Two  dimensional  row  and  column  vectors 


Subroutine:  ACTU 


Description:   Composes  the  actuator  matricies  (Eqs.  3.42  and 
3.43) 


Calling  arguements:   CALL  ACTU( Fa, Ga , Ha , Nf a , Nga , Nha ) 

Input  arguments:   None 

Output  arguments: 

Fa,  Ga,  Ha  Actuator  matricies 

Nfa,  Nga,  Nha       Two  dimensional  Row  and  column  vectors 
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Subroutine:   PLANT 


Description:   Composes  the  basic  airframe  plus  actuator 
matricies  (Eqs.  3.44  and  3.45) 


Calling  arguments:   CALL  PLANT ( Fm, Gm, Hm, Dm, Fa , Ga , Ha , Fp, Gp, Hp, 

Nfm, Ngm, Nhm, Ndm, Nfa, Nga,Nha, Nfp, Ngp, Nhp) 


Input  arguments 


Fm,  Gm,  Hm,  Dm         Modified  basic  airframe  matricies 

(Eqs.  3.44  and  3.45)   Generated  in 
subroutine  MODEQ. 

Fa,  Ga ,  Ha  Actuator  matricies  (Ecs.  3.42  and  3.43) 

Generated  in  subroutine  ACTU. 

Nfm,  Ngm,  Nhm,  Ndm     Two  dimensional  row  and  column  vectors 
Nfa,  Nga,  Nha 


Output  arguements: 

Fp,  Gp,  Hp  Airframe  plus  actuator  matricies 

Nfp,  Ngp,  Nhp  Row  and  column  vectors 
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Subroutine:   LAVS 


Description:   Computes  the  function  gains  (App.  A) 

Computes  the  coefficients  in  the  control  law 
matr icies 

Composes  the  control  law  matricies 


Calling  arguments:   CALL  LAWS ( alpha, nz , psi , qc, ri , ts , Ac , Bfc, 

Be,  Cc, Dfc,Dc,Nac,Nbfc,Nbc,  Ncc,Ndfc,  Ndc) 


Input  arguments: 

alpha  Steady  state  angle  of  attack  in  degrees 

nz  Steady  state  normal  acceleration  =  1.0 

psi  standard  atmosphere  pressure  at  steady  state 

altitude  in  lb/ft**2 

qc  dynamic  pressure  in  lb/ft**2 

ri  pressure  ratio 

ts  sampling  time  in  seconds 

Output  arguments: 

Ac,  Bfc,  Be       Control  law  matricies  (Eqs.  3.29  and  3.30) 
Cc,  Dfc,  Dc 

Nac,  Nbfc,  Nbc    Row  and  column  vectors 
Ncc,  Ndfc,  Ndc 
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Subroutine:   VC7AIN 


Description:   Computes  the  variable  gain  matrix  based  on 
control  surface  impairment.   The  subroutine 
is  disigned  to  use  reconfiguration  algorithm 
introduced  in  Ref.  1. 


Calling  arguments:   CALL  VGAINf i fai 1 , i fix , Gm, GmO , GAIN, Nam, 

NgmO,  Ngain) 


Input  arguments: 

ifail         Control  surface  failure  flag  set  in  subroutine 
•LONLAT;       No  failure  =  0 
Failure  =  1 

ifix  Fix  parameter.   The  GAIN  matrix  will  be  computed 

fo  the  impaired  aircraft  if  ifix  =  1  and 
ifail  =  1.   Ifix  is  read  interactively  by  sub- 
routine FLITE1 . 

Gm  Airframe  control  matrix.  From  subroutine  MODEO. 

GmO  Unimpaired  airframe  control  matrix.   Generated 

in  subroutine  LONLAT.   Note  if  none  of  the 
control  surfaces  are  damaged  or  failed,  Gm=GmO . 

Ngm,  NgmO     Row  and  column  vectors 


Output  Arguements: 

GAIN         Variable  gain  matrix  described  in  section 
III. A. 

Ngain         Row  and  column  vector 
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Subroutine:  EXPINT   [Ref.5] 

Description:   Computes  both  the  matrix  exponential 

Fps*ts 
Aps  =   e 

and  the  integral 

>  Fps*s 
\    e      ds 

Reference  5  gives  the  method  used  to  perform 
the  computations. 


Calling  arguements:   CALL  EXPINTf Fps , Nfps , Aps , Naps , Duml , 

Nduml  ,  ts  ,  iop,  Dum.2  ) 


Input  arguments: 

Fps  Airframe  plus  actuator  plus  sensor  transfer 

matrix.   Generated  in  ARCRFT 

Nfps         Two  dimensional  row  and  column  vector 

ts  Samplina  time 

iop  Print  parameter:  0  =  Do  not  print  results 

Fp s  *  t  s 
Otherwise  print  Fps,  Ts,  e 

Dum2         Vector  of  working  space  for  computations 

Output  arguments: 

Aps         Matrix  exponential  described  above 

Naps        Two  dimensional  row  and  column  vector 

Duml       'The  integral  of  the  matrix  exponential.  To  be  used 
in  computing  the  discrete  control  matrix  Bps. 

43x10    43x43    43x10 
Bps   =   Duml  x   Gps 
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The  following  subroutines  perform  the  matrix  operations 
necessary  to  compose  the  A/F-18  system  matricies. 


Subroutine:  MULT  [Ref.5] 

Description:   Performs  the  matrix  multiplications. 

Calling  arguments:   CALL  MULT( A, NA, B, NB, C, NO 

Input  arguments: 

A,  B         Matricies  to  be  multiplied 

Na ,  Nb       Two  dimensional  vector  giving  the  number  of  rows 
and  columns  in  A,  and  B. 

Output  arguments: 

C  Product  of  A  and  B,   C  =  AB 

Nc  Row  and  column  vector 

Subroutine  ADD  and  SUBT  are  called  just  as  MULT. 
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Subroutine:   NULL   [Ref.  5] 
Description:   Generates  a  null  matrix 
Calling  arguments:   CALL  NULL(A,NA) 

Input  arguments: 

Na  Two  dimensional  vevtor  givinq  the  number  of  rows 

and  columns  of  the  null  matrix, 

Na ( 1 )  =  Number  of  rows 
Na ( 2 )  =  Number  of  columns 

Output  arguments: 

A  Null  matrix  havinq  dimension  Na ( 1 )  x  Na ( 2 ) 

Subroutine:   JUXTR  [Ref.  5] 

Description:   Juxtaposes  two  matricies  by  row 

Calling  arguements:   CALL  JUXTR( A, Na , B, Nb, C, Nc ) 

Input  arguements: 

A,  B        The  two  matricies  to  be  juxtaposed  by  row 

Na ,  Nb      Row  and  column  vectors 

Output  arguments: 

A 
C  The  juxtaposition  of  A  and  B    C  =  

B 
Nc  Row  and  column  vector 
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Subroutine:   JUXTC   [Ref.5] 

Description:   Juxtaposes  two  matricies  by  column 

Calling  arguments:   CALL  JUXTR( A, Na , B, Nb, C, Nc ) 

Input  arguements: 

A,  B        The  two  matricies  to  be  juxtaposed  by  column 

Na,  Nb      Row  and  column  vectors 

Output  arguments: 

C  The  juxtaposition  of  A  and  B    C  =  A  I  B 

Nc  Row  and  column  vector 


Subroutines:   ADD  and  MULT 

The  calling  seguences  for  these  subroutines  are  exactly  as 

described  for  those  matrix  operations  above. 
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Subroutine:   OUTPUT 

Description:   Outputs  a  matrix  with  rows  and  columns 
numbered . 

Calling  Arguments:   CALL  OUTPUT ( A, Na f 1 ) , Na ( 2 ) , ' A    ') 

Input  arguments: 

A  Matrix  to  be  output 

Na  ( 1 ) ,  Na ( 2 )         Number  of  rows  and  number  of  columns 

in  A  respectfully 

'A  Matrix  name.  Must  be  placed  in  guotes 

and  occupy  four  spaces. 

Output  arguments: 

The  A  matrix  is  written  to  the  FA18  RESULTS  file 
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APPENDIX  I. 
F/A-18  EXEC  PROGRAM 


&TRACE  OFF 

CLRSCRN 

GLOBAL  TXTLIB  VALTLIB  VFORTLIB  IMSLDP   NONIMSL   CMSLIB 

GLOBAL  LOADLIB  VFLODLIB 

FILEDEF  06  TERMINAL 

FILEDEF  01  DISK  F/A-18  DATA  Al 

FILEDEF  02  DISK  F/A-18  RESULTS  ( LRECL  132 

FILEDEF  03  DISK  OPTMATD  DATA  Al 

FILEDEF  04  DISK  OPTPLOT  DATA  Al 

&TYPE  LOADING  FA18  SIMULATION 

LOAD  FA18  ORACLS1  (START 

&EXIT 
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APPENDIX  J 
FLIGHT  CONDITIONS  AND  STABILITY  AND  CONTROL  DERIVATIVES 


FLIGHT  CONDITIONS 

MACH  =   0.6 
ALPHA  =  2.6184 
PCTHLC  =  20.496 
DLHTD  =  0.68914 
DLAD  =0.0 
DLEFL  =  3.0968 
DTEFL  =  3.2647 
DRUDL  =0.0 
FLARUD  =0.0 
CG  =  0.2206 
XIYY=  0.1220E+06 


ALT  =  10000.0 

VTRUE  = 

=  646.42 

THRUST  =  3531.5 

GAMA  = 

0.0 

PCTHRC  =  20.496 

XIXX  - 

21227.0 

DRHDT  =  0.68914 

XIXZ  = 

-1827.0 

DRAD  =0.0 

DLEFR  =  3.0968 

DTEFR  =  3.2647 

DRUDR  =0.0 

DELSB  =0.0 

WAIT  =  32550 

XIZZ  =  0.13958E+06 

B.   LONGITUDINAL  DERIVATIVES 


XU  =  -0.13257E-01 

XW  =  0.71265E-01 

XQ  =  0.32650 

XWD  =  0.39729E-03 

XDSB  =0.0 

XDTH  =  0.14257 

XDS  =  -0.43034E-03 


ZU  =  -0.73337E-01 

ZW  =  -1 .1526 

ZO  =  -5.6525 

ZWD  =  -0.6917 

ZDSB  =  0.0 

ZDTH  =0.0 

ZDS  =  -0.19801E-01 


XDLF  =  -0.412942E-03   ZDLF  =  0.64114E-02 


MU  =  -0.12988E-04 

MW  =  -0.11331E-01 

MO  =  -0.59346 

MWD  =  -0.34049E-03 

MDSB  =0.0 

MDTH  =0.0 

MDS  =  -0.24151E-02 

MDLF  =  -0. 5230E-03 


XDTF  =  0.17674E-03 


ZDTF  =  -0.45250E-01   MDTF  =  0.32773E-01 


C.   LATERAL-DIRECTIONAL  DERIVATIVES 


YV  =  -0.2437 

YVD  =0.0 

YR  =  0.77567 

YP  =  0.67478E-02 

YDA  =  0.86683E-03 

YDR  =  0.52918E-02 

YDLF  =0.0 

YDHT  =  -0.21027E-02 

YDTF  =0.0 


LV  =  -0.29125E-01 

LVD  =0.0 

LR  =  0.70076 

LP  =  -3.1241 

LDA  =  -0.43808E-02 

LDR  =  0.72058E-03 

LDLF  =0.0 

LDHT  =  0.39478E-02 

LDTF  =  0.24948 


NV  =  0.69481E-02 

NVD  =0.0 

NR  =  -0.21705 

NP  =  -0.15254E-01 

NDA  =  0.77305E-04 

NDR  =  -0.46503E-03 

NDLF  =0.0 

NDHT  =  0.57361E-04 

NDTF  =  -0.32655E-02 
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FLIGHT 


DEFINITION 


UNITS 


XMACH 

ALT 

VTRUE 

ALFA 

THRUST 

GAMA 

PCTHLC 

PCTHRC 

DLHTD 

DRHTD 

DLAD 

DRAD 

DLEFL 

DLEFR 

DTEFL 

DTEFR 

DRUDL 

DRUDR 

FLARUD 

DELSB 

CG 

WAIT 

XIXX 

XIYY 

XIZZ 

XIXZ 


MACH  NUMBER  ND 

ALTITUDE.  MSL  FT 

TRUE  AIRSPEED  FT/SEC 

ANGLE  OF  ATTACK  DEG 

TOTAL  ENGINE  THRUST  LB 

FLIGHT  PATH  ANGLE  RAD 
LEFT  THROTTLE  POSITION  COMMAND  (100K-MAX    AB)         PERCENT 
RIGHT  THROTTLE  POSITION  COMMAND! I OOX-MAX    AB)        PERCENT 

LEFT  HORIZONTAL  STAB.  DEFLECTION  (+  TE  DOWN)  DEG 

RIGHT  HORIZONTAL  STAB.  DEFLECTION  (+TE  DOWN)  DEG 

LEFT  AILERON  DEFLECTION  (+TE  DOWTJ)  DEG 

RIGHT  AILERON  DEFLECTION  (+TE  DOWN)  DEG 

LEFT  LEADING  EDGE  FLAP  DEFL  (HE  DOWN)  DEG 

RIGHT  LEADING  EDGE  FLAP  DEFL.  (+LE  DOWN)  DEG 

LEFT  TRAILING  EDGE  FLAP  DEFL.  (+TE  DOWN)  DEG 

RIGHT  TRAILING  EDGE  FLAP  DEFL  (+TE  DOWN)  DEG 

LEFT  RUDDER  DEFLECTION  (  +  TE  LEFT)  DEG 

RIGHT  RUDDER  DEFLECTION  (+TE  LEFT)  DEG 

FLARED  RUDDER   (  PA/TO  ONLY  )  DEG 

SPEEDBRAKE  DEFLECTION  DEG 

LONG.  CENTER  OF  GRAVITY  LOCATION  IN  PCT  MAC  PERCENT 

AIRCRAFT  WEIGHT  LB 

MOMENT  Cf  INERTIA  ABOUT  X-BODY  AXIS  SLUG-FT2 

MOMENT  OF  INERTIA  ABOUT  Y--BODY  AXIS  SLUG-FT2 

MOMENT  OF  INERTIA  ABOUT  Z-BOOY  AXIS  SLUG-FT2 

XZ  PLANE  PRODUCT  OF  INERTIA  SLUG-FT2 
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APPENDIX  K 
F/A-18  RESULTS  FILE 
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